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1. INTRODUCTION 


1.1 Objective : The objective of this work is to develop capabilities to 
evaluate the Wiener Spectrum and Modulation Transfer Function of the 
autoradiographic enhancement process. The autoradiographic enhancement 
process is a technique for increasing the amount of information which is 
recovered from photographic film, plates and prints. The technique consists 
of radioactivating the image silver by chemically combining ic with a beta 
emitting isotope, sulfur-35; the autoradiograph, an enhanced copy of the un- 
derexposed imag'?. is made by exposing a second film to the radiation from 
the underexposed anginal. The efficient response of the photographic 
emulsion to beta radiation produces an autoradiograph image which is enhanced in 
density and contrast so that images which were invisible or only faintly visible 
on the original film can be easily seen on the autoradiograph. 

The Perkin-Elmer 1Q10A Microd*4nsitometer is used to generate raw data 
from film. 

Evaluation of the Wiener Spectrum and Modulation Transfer Function provides 
a measure of effectiveness of image recovery in terms of signal-to-nolse ratio 
and detective quantum efficiency. 'It is the quantitative evaluation of the 
information transfer due to autoradiographic intensification. 

Software is being developed to provide data analysis capabilities for 
autoradiographic image processing and other image processing applications where 
the microdensitometer is being used. The data analysis capabilities shall be 
developed on the NOVA 2 compu ter located in room 369 of building 4481. While 
IDAPS provides many of the analysis capabilities to be developed on the NOVA, 
IDAPS is designed for only 8-bits of resolution while the microdensitometer output 
provides 12-bits of resolution. For the autoradiographic enhancement process, the 
12-bits of resolution Is definitely needed, especially for signal to noise ratio 
calculations. Also, the data on microdensitometer output tape cannot be used as 
is with IDAPS because of header .information included in each data record and 
because of the 12-bit length of the data. 

The effort reported in this document is the software requirements definition and 
the preliminary software design of the analysis capabilities to be implemented on 
NOVA 2. 
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2. SYSTEM CONFIGURATION 


Figure 2.1 illustrates tlie system configuration for data analysis. 
T.ne piece of film is scanned by the microdensitometer and the image density 
is recorded on a 9-tracK, 800 BPI magnetic tape. Tne tape is input to the 
NOVA 2 w*"ire the data is analyzed and a copy of the results Is obtained from 
the Textronix 4014 hard copy unit. 
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i. rgp S QKI iV.-VJi 


The software shall be designed to analyze the film-scan data obtained 
from the Perkin- Elmars Mi crodensitometer. 

1*1 Inputs 

The Inputs to the software shall be: 

(1) The microdensitometer output tape 

(2) User-specified variables 

_LuU L.MlCltfdr nanometer Ou t put Tape 

The icia from the mlcrodensltom^ter output tape Is obtained via a 9-trac’c 
magnetic tape. The taoe da*a format Is Ulustra*ed in Figure 3.1. 

A nugnetic tape record Is generated by the .nicrodensitometer control 
system at the conclusion of each scan line, or segment thereof. Tsrminatior. of 
the record consists of an IBM -compatible Inte^-Record-Gap. 

after the last record has been written for a given scanning program, the 
mlcrodensltomc-ter control system causes the writing of an IBM -compatible standard 

End of File Mar*. 

lAHS F<?lf!K>t; 

1. Start o f Record - Two tape charact*rs containing the hexidecimai characters 

"F t FT" . These characters have no relation or value relative to any of the data to 
follow. 

2. Identification - Forty tape characters, each containing an eight-bit byte representing 
the ASCII code for the alphanumeric characters entered during the IDENT mode of 
operation. 

3. Numbe r o f Samples - Two tape characters containing a sixteen-bit binary 
number representing the (positive) number of data values included in the record, 
following completion of the coordinate information. 
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T aua format (cont'd)- 


4. X-Coordina ie - Four tape characters containing a thl'ty-two bit Binary number 
representing the position ulong the X-axls at which the current scan line commenced. 
The number is justified to the LSB end of the word; a negative position is Indicated 
as a two's complement of the absolute value. All position readings are In microns. 

5. Y- Coordinate - "our tape characters containing a thirty-two bit binary number 
representing the Y-axls coordinate at which the current scan line commenced. 
Interpretation as noted above. 

6. Delta X - Four tape characters containing a thirty-two bit binary number 
representing the sampling Interval along the X-axls. Interpretation as above. (A 
negative number Indicates that the direction cf scanning was negative.) 

7. Data Val ues in Buf fer - Two tape characters for each of the data values 
enumerated in "3" abo/e. Each pair of characters contains a sixteen-bit binary 
number, the leading twelve bits of which are significant. The final four bits are 

to be discarded. Interpretation of the data Is as follows*. Convert the binary numLer 
to decimal and divide by 800. The quotient is the data value for the density or 
transmittance reading, as recorded dur .g the scan. 

The data corres pending tv, a scan line is stored on one record provided 
the number of samples is no greater than 3200. If the number of samples per scan 
line is greater than 3200, then the scan line data Is segmented and stored on more 
than one record on tape. Start of each new i can line, however, causes a new 
record tc be written. 


3.1.2 Use r-Sieclfl ed Variables 

T vj following user-specified variables are supplied In response tc 

software-query: 

(a) No. of File NF, to be read from the microdensltornoter tape. 

(b) The Diffuse Density Factor, DD 
vc) The Maximum Lag Factor, M 

(d) The Length of the scanning slit, SL 


I-f. 


U_AualYam vimaiimqi 


Tie software shall be designed to have the following da? a analysis 
capabilities: 

(1) Compute Histograms 

U) FllUr Data 

(3) Compute Means 

(•1) Coinjjte Sta.oda-d Deviation 

(5) Compute Autocorrelation functions 

(o) Compute Power Spectral Dsi ilty Functions 

The software shall be designed to permit future capabilities to bo 
developed and integrated with the existing software. 

3.2.1 H istograms 

The density data range from a scan frame shall be divided into fifty 
classes. The data from a scan frame shall be distributed ovor these classes to 
generate a histogram. 

Filters. 

Tns software shall have the capability to filter raw data from the 
nlorodensltometer. The specific filter or filters are yet to be determined. 

3.2.3 M^ans 

The software shall have the capability to compute mean of the raw data 
and mean of the data normalized by a diffuse density facto:. The software shall 
normalize data to zero mean. 

3. 2.-1 Standard Deviation 

The software shall have the capability to compute mean square value and 
standard deviation of the data 
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2^.5 A i»ocom ? l fl » i i 0 p Function? 

Using the user-specified maximum Is? factor, the software shall be 

3»in.° ^r putln9 . a . ,,OCOrre i atl0n luncMons a ‘*o of normalizing ajtocortelauons 
unctions with respect to zero-lag autocorrelation function. 

Power Soectral E >?n« lty 

It shall aUn ™ tW ! r * 3haU b ® MfVi;3le of computing power spectral density functions, 
it shell also compute smoothed power spectral density functions. 

1.3 Outputs 

The software shall have the following outputs: 

(1) Data Identification from Microdensitometer Tape 

(2) Raw data mean 

(3) Mean Square Value 

(4) Standard Deviation Value 

(5) Histogram table 

(6) Autocorrelation and Normalized Autocorrelation Functions 

(7/ Power Spectral Dsr.sity and Smoothed Power Spectral Density Functions. 
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4. PRELIMINARY SOFTWARE DESIGN 


Tho software design 1$ modular. An overview of the preliminary software 
design showing the top-modules Is illustrated in Figure 4.1. The main program calls 
Independent modules which Implement unique data in; at, data preparation and data 
analysis capabilities. The sequence of execution Is top-down and left to right. 

4U— Spurting. Quarts 

The detailed software design is presented in structure charts shown in 
Figures 4.2 thru A. a and the HIPO's of section 4.2. T ie sequence of execution 
of modules in the structure chart Is ton-donn and lc t to .4ght. Each module has one 
entry and one exit. Control always returns ‘o the next statement In the calling 
module. Arrows on a line indicate multiple calls. A»tows on a group of lines indicate 
multiple sequence of calls. 
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Figure 4.2 Partial Structure Chart Data Preparation Modules 




































HIPO'S 

The Hierarchical Input-Process-Cutput (HIPC) charts are given for the 
modules shown on the structure charts. A Hl a O gives an overview of the actions 
performeJ vtthln a -nodule. 


4.2.1 HIPO-MAiN 


PURPOSE: Mi in accepts the user-supplied data and calls other modules for 
preparation of tape data and data analysis 


Input 


User 

supplied 

data 

NF, tape 


RD, N? , 
A, B, NC 


Process 


Ask for and Read from programmer's console the 
File No. NF, the diffuse density factor KDD, the 
max lag factor M and length of scanning slit L. 

Call (a module) PTDT to prepare Jljl: file from tape 
data and to get H ■ Delta X. 

Call HIST to prepare histogram 


NF, KDD 
M, L 


RD, NP, 
H 


RD, NP Call F1LT to filter raw data 

RD, NP , Call MEANS to compute mean, convert data to density 
KDD units and to normalize data to zero mean 

ZRD, UP Call SIGMA *.o compute standard deviation and Mean 
Square Value 

ZRD NP , Call A'. TOCOR to compute (M + 1) autocorrelation 
M functions and to noimalize w.r.t. autocorrelation 

with r = 0. 


RRH, H, L 
M 


Call PS Df to compute 5^ AND S k , 


GKC, 

GKH 


STOP 










l.J .2 HI F ?■ PTDT 


Purpose: PTDT prepares disk (lie o( raw data from tape data. 


Step No, 

Input 

Process 

Output 

1. 

Ta pe 

Call MTOPD to position tape at device MT0 


*■» 

NF 

If (NF.GT. 1) call SKPFL to skip files on tape to come 
to file no. NF 


3. 


Call RDFL to read File No. NF from tape Into disk 
file TD. NR Is number of records. NW ts array of 
no. of words in each record. 

TD, NW, 
NR 

4. 


Call HDR to print header information and get delta 
X - NH 

NH 

5. 

TD, NW 
NR 

Call RFMT to reformat data lr. disk File TD to raw 
da.;* in disk file RD 

RD. NP 

6 • 


RETURN 
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Purports: SiOTL skips Tiles on tape to ic-a.l tile no. NT. 


Input 


Process 


1. NF 


Call MTDIO to read one record from tape and word- 
count ICNT in that record. 


ICNT 


If ICNT is greater than 1 , p to step i 


iCNT 


If ICNT « 1, this indicates EOF. 


Increment file counter KFL. 


If KFL » NF, RETURN 


Else go to step 1 . 
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-1.2.4 HiPO-RDFL 



Step No. Input 



ICNT 


INPBUF 


Process 


Call FOPEN to open a disk file TD to receive tape 
data 

Call MTDiO to read one record from tape into 
Temporary array INPBUF 

If ICN7 * i (EOF), then go to step 6. 

Else record number of records, no. of words per 
rscord and call WRITR to put data in file TD. 


Output 


INPBUF, 

ICNT 


Go to step 2 


Return 
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4.2.5 H1P0-HD3 


Purpose: HD4 r»ads llrst record In dis'<. (lie TD and prints I.D. and computes 
H - Dt ita X. 



■sn 

Process 

Output 

1 . 

TD 

Ceil READR to read first record rom IlsX file TD 
to Temp, array ITEMP 

ITEMP 

2. 

ITEM? 

Print ID *.n words 2 through 21 of ITEMP 


3. 

ITEMP 

Print no. of samples * HEMP (22) 


4. 

ITEMf 

Call COPD 3 times to compute X, Y coordinates 
and D;ita X. After each call, print the coordinates 
and Delta X - NH. 

NH 

5. 


Return 
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f — ■ 


Purpose: CORD computes X. Y coordinates onu Delta X from 2 input words 

for each value 


Step NoJ Input 


Process 


ITEMT, I Load a double precision word IWD with 2 input 
words 


IWD 


If an input word Is negative, use specia’ techniques 
'or loading Into IWD. 


Return 


-2 




4.2.7 H1P0-RFMT 

Purpose: RTMT reformat! samp)* values lata in disk file TD and loals the 
data values only into disk file RD. 


Process 


Output 



ITEM- , 
NR, NW 


JTEMP, 

NF* 


Call FOPEN to open disk file RD 

Call READ4 to read a record from file TD into 
.emporary array ITEMP 

Call GBYTE io get 12 most significant bits from 
each dn*a point (words 29 tnrough NW 'or each of 
NR records) and store them 4 .n Temp, array JTEMP. 

Call WR»TR to write JTEMP in file RD. 


Repeat steps 3 and 4 NR times 

Compute No. of data points NP In each record 

READ j record from file RDand print it for 
verification. 


ITEMP 


JTEM? 


ITEM? 


RETURN 
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4. ’.6 HIPO-HI37 


Purpose: HIST computes the histogram for the raw data. 


Step No. 

Input 

Process 

Output 

1. 

RD, NW 

Call READ* to read a record from RD file Into 
temporary array ITEMP 

ITEMP 

2. 

ITEMP , 

ic, np: 

Call DClS to distribute data in ITEMP into class 
Intervals. 

NHST 

3. 

NHST 

Sum NHST into MST by repeating steps 1 and 2 for 
NR records 

MST 

4. 

M3T 

Find IMIN and IM\X such that MST has nonzero 
entires between IMIN and I MAI; 

IMIN, 

’.MAX 

5. 

IMIN , 
IMAX, 
M3T 

Type MST values 


6. 


Return 
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9 hpo-dclS 


Purpose: EXJi.S assigns eash data value to a class Interval fcr calculating 
histogram. 


Step No, 


1 . 


Input 


ITEV1P, 
IC. N PI 

I 


Process 


Compute Integer ^art of (ITEM? (J)/l C) * I for one 
value of J. 

Increment NH3T (i) by 1. 

Repeat steps 1 and 2 for NPI points 
Return 
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4.2.10 HIPO-MEaNS 


Purpose: MEANS calls other modules to compute mean of raw data and normalize 
data to diffuse density units and ^ero mean. 


Step No. 

Input 

Process 

1 

Output 

1. 

ICH, NR, 
NP 

Call AVG to compute raw mean. 

RA\VM 

2. 

KDD 

Compute the diffuse density multiplier DDM '-om 
KDD as DDM = KDD/RAWM 

DDM 

s » 

DD:Vj , 

ICHK, NR, 
NP, ICHN 

Call NRM1Z to normalize data to diffuse density 
units . 


4. 

KDD 

F r int mean of normalized data KDD. 


5. 

ICHN, 
ICHZ, NR, 
NP, KDD 

Call ZMEAN to normalize data to zero mean 


6. 


Return 
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4.2.11 H1PO-AVG 


Pu.pose: AVG commutes the mean (average) of a group of records on a disk file. 



warn 

Process 

Output 

1. 


Call READP. to read a record from file RD into Temp, 
array ITEMP. 

ITEM? 

2. 

NP, ITEMP 

Compute cumulative sum IS JM of points in ITEM? 

ISU M 

3. 

NP 

Compute cumulative sum of number of points in 
aach record. 

INP 

4. 

NH 

Repeat steps 1 , 2 and 3 for NR records 

I SUM 

5. 

INP, ISU N 

Compute M.- ; an RAWM * ISU M /INP 

RAWM 

6. 

RAWM 

Print RA VM 


n 

• • 


Return 
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•4.2.12 HIPO-NRMLZ 


Purpose: NRML2 normalizes data In a disk file by multiplying It with a factor 
and stores it in another disk file. 


Step No.{ 

1 . 

2 . 

3. 

4 . 

5 . 

6 . 


Input 

ICHN 

DDM, 

TEMP 

NR 


Process 

Gall FOf-EN to open a file on disk 
Call READ< to read a record Into Temp array 
Multiply each of NPI data values by DDM 
Call WrtITR tr. write TEMP in disk file NRD 
Repeat steps 2 , 3 and 4 for NR '■ecords 
Return 


Output 


TEMP 

TKM? 
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4.2.13 H1PO-ZMLAN 


Purpose: ZMEAN normalizes data on a disk file by subtracting the data mean 
factor from each value. It s'ores th® resulting data in another disk 



Process 


1. ICHZ 


Jail t'OPEN to open a file on disk 


2. ICHN 


Call READR to read a record into temporary array 
TEM? 


TEMP 


KDD, 

W?I 


Subtract KDD from each of NPI data values 


TEMP 


TEMP 


Call WRITR to write TEMP into a disk file ZRD 


Repeat steps 2, 3 and 4 for NR records 


Return 





-L2.N HIPO-SIGMm 

Pjrpose: SIGMA .imputes the mean square value and root mean square value of 
data In a disk file. 


Step No. 

Input 

Process 

Output 

1. 

ICH 

Call READR to read a record from file ZRO Into 
Temp array TEM? 

TEMP 

*• • 

NPI, 

Compute cumulative sum of squares, SSQ for NPI 

SSQ 


TEM? 

data values In TEMP 

3. 

NP 

Compute sum of no. of values for NR ecords 

INP 

4. 

NR 

Repeat steps 1 , 2 and 3 for NR records 

SSQ 

5. 

INP, ssg 

Compute XSQ * 33Q/INP 

XSQ 

C. 

XSQ, INP 

Compute SIG * (XSQ*NP/ (NP-1) ) 1/2 

SIG 

7. 

XSQ, SIG 

Print XSQ and SIG 


8. 


Return 
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4.2.15 H1PO-AUTOCUR 

Purpose: AUTOCOR computes the autocorrelation and normalized autocorrelation 


Step No. 

Input 

Process 


1. 


Set XT, RAT, TFMP to zero 

XT 

2. 

ICHZ 

Call READP> to read a -eiord from file ZRD Into Temp 

TEMP 

3. 

NP 

Compute cumulative sum INP 

INP 

4. 

TEMP, M 

Shift data value contents of TEMP down by M 

TEMP 

3. 

M, XT, 
TEMP 

Transfer contents of XT Into first M values of TEMP 

TEMP 

6. 

M, XT, 
TEMP 

Transfer last M data values In TEMP to XT 

XT 

7 . 

RAT, TEMl 

1 For JR * 0 through M, compute cumulative sum for 
I - 1 to NPI,RAT (JR+1) - RAT (JR+1) + TEMP (I) * 
TEMP (I + JR) 

RAT 

8. 

NR 

Repeat steps 2 through 7 for NR records 


9. 


Set TEMP to zero and do steps 5 and 7 one time 

RAT 

10. 

RAT, MR 

Compute RAT 'J)/(INP-J+1) for J = 1 , MR 

RAT 

11. 

RAT, MR 

Compute normalized autocorrelations function 
RATN (J) = RAT (J)/RAT (1) for J * 1 through MR 

RATN 

12. 


Print J, RAT(J) and RATN (J) 


13. 


Return 
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4 . 2 16HIP0 : LFF 


Purpose : 


LFF is the low frequency filter applied to raw data 
f andard deviation, autocorrelation and power 
spectral density functions are computed. 


Step No 

Input 

Process 

HiirmiP 

1 

NP.P 

Compute number of records in which the 
first (2P+1) points lie. 

NREC 

2 

P.NP 

Compute the record number in which the 
(P+l)th point lies. 

QUO 

3 

TEMPI , P 

Calculate sum of first (2P+1) points. 

SUM 

4 

QUO , P , NP 

Calculate location of (P+l)th point in 
record no. QUO. 

LOCI 

5 

TEMPI, 
LOCI, SUM 
P 

(P+1) th point as TEMPI (LOCI) - 
TEMPI (LOCI) - SUM/(2P+1) and store TEMPI 
on disk. 

TEMPI 

6 

NP.N.P 

For I - (P+2) to (N-P), where N is the 

total number of points, compute record 
nos. in which points I, I-P-l and I+P lie. 

RECI, 

REC2, 

REC3 

7 

REC1 , 
REC2, 
REC3, NP 

Calculate location of points I, I-P-l and 
I+P in their respective records. 

LOCI, 

L0C2, 

L0C3 

8 

reci, rec; 

REC3 

.Read respective records from disk. 

TEMPI , 
TEMP 2 , 

9 

LOC2.LOC3 
TEMP 2 , 
TEMP 3 

Calculate SUM - SUM-TEMP2 (L0C2) + TEMP3(LOC3; 

TEMP3 

SUM 

10 

TEMPI, LOCI 
SUM, P 

Calculate XT - TEMPI (LOCI) - SUM/(2P-t-l) 
and store on disk at TEMPI (LOCI). 

XT, TEMPI 

11 


Repeat steps 6 through 10 till done. 


12 


Return . 

1 

1 
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^•2.17 HIPQ.HF T 

Purpose : HFF is the high frequency filter applied to output 

data of the low frequency filter. 


Input 


Process 


NP.Q 


Q.NP 


QUO.Q,, 

NP 

Q , SUM 


NRN.Q 


REC1 , 
REC2 , 
REC3.NP 

REC1, 

REC2, 

REC3 

L0C2 , 
L0C3 , 
TEMP2 , 
TEMP 3 

SUM , Q 


Calculate number of records in which the 
first (2Q+1) points lie. 

Compute the record numb ’' in which the 
(Q+l)th point lies. 


NREC 


TEMPI ,Q Calculate sum of first (2Q+1) points 


Calculate location of (Q+i)th point in 
record No. QUO. 

Calculate new (Q+l)th point as 
TEMPI (LOCI) - SUM/ (2Q+1) and store 
TEMPI on disk. 

For I ■ (Q+2) to (N-Q) . where N is the 
total number of points, compute record 
numbers in which points I, I-Q-l and 
I+Q lie. 

Calculate location of points I, I-Q-l and 
I+Q in their respective records. 


Read REC1, REC2 and REC3 from disk. 


Calculate SUM - SUM - TEMP2(LOC2) + 
TEMP 3 (L0C3) 


Calcuiatt XT * SUM/(2Q+l) and store on 
disk at TEMPI (LOCI) 


LOCI 

TEMPI 

(LOCI) 


RECl . 
REC2, 
REC3 


LOCI, 
L0C2 . 
L0C3 

TEMPI, 
TEMP 2 , 
TEMP 3 


XT , TEMP 


Repeat steps 6 though 10 till done 


12 


Return . 












4.2.19 HIPOINF 


Purpose INF is the Inverse filter used on power spertral 

density function as computed from filtered data. 



wmm 

Process 

Output 

1 

GKC, Q 

For K - 0, Set , 

GKCP (K+l) - <2Q+1) Z 

GKCP(l) 

2 

GKC.Q.M 

For K - 1,2 M, Compute 




GKCP (K+l) - GKC(K+i) • 

r ((2Q+l) • Sin(Kir/2M)) / 
Sin ((2Q+J.) K*/2M)/] 2 

GKCP 

3 


Return. 
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*.2.20 HIPO AC 

Purpose AC is the Aperture Compensation applied to the power 
spectral density data output from INF. 


Process 


GKCP 


GKCP.M.H, 

XA 


For K-0, Set 
GKC3P(K+l) - 0 


For K-1,2 M, compute 

GKC3P(k*l) • GKCP(K+1) • 

(k-tt • XA/2H- M* Sin (k • n • XA/2H- M) ) 


GKC3P(1) 


GKC3P 


Return. 






I J 




l HIPO - SPDF 

Purpose : SPDF compute, the smoothed power spectral density function, 

Fiom the power spectral density functions. 


Input 


Process 


GKC.MR 


KR.GKC, 

GKH 


Compute 

GKH(l) - 0.5 t GKC(i) + GKC(2)J GKHi 

GKH (MR) - 0. 5 fGKC(MR-l) + GKC(MR)J GKHi 

For KR»2 (MR -1), compute 

GKH(KR) - 0. 25(CKC(KR-1) ♦ GKC(KR+1)) + 0.5 CKH 
GKC(KR) 

Print KR, GKC(KR), GKH(KR) 


Return 







APPENDIX A 


Appendix A contains the mathematical formulae and notations used 
for analysis of autoradiographic image data. These mathematical 
formulae form the basis for implementation into programs. 

Notation 

1 ■ Slit height 

a - Slit width 

A or H - Sampling Interval 

D • Raw Data Point 
n 

N - Total number of data points 

B * Mean of raw data, D„ 

n n 

SIG ■ Standard Deviation 

d - Raw Data converted to diffuse density 
n 

Z n - Zero Mean Data (in diffuse density) 

Z^ - Zero Mean Data after low frequency filtering 

Z ^ ■ Zero Mean Data after high frequency filtering 

p * Parameter in low frequency filter 

q - Parameter in high frequency filter 

R(r) - Autocorrelation function 

r - lag factor 

m * Maximum lag factor 

DDF * Diffuse density factor 

* Spectral density function 

G£. ■ Spectral density function after inverse filter 

■ Spectral density function after aperture compensation 

* Smoothed spectral density function 

k «* 0.1,2....m 

f ■ Frequency « k 

2"S~m 


User-specified variables. 
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Raw Mean 


5 «'S n i °n 


Normalized Data 


U - D / DDF 
n n 


Zero-Mean Data 


Z - U - tJ 
n n n 


Mean Square Value 
MSV - 


* " d 2 


Standard Deviation 


SIG - 


f n ! 

( 5PT MSV ) * 


Low-frequency Filter 


Z ' • Z 

4 n 4 


n * 25*1“ j£. p - 2 J 


+n 


High Frequency Filter 
Z 


n 


1 

7q+r 


q 

z 

J— q 


Z' . 
J+n 


Autocorrelation Function 


K 


1 N-r 
iPr Z 
n-1 


'n * Z 


n+r 


n « p+1 ,p+2 , . . ,N-p 
p * appx. 100 to 1000 

n - q+l ,q+2 .... ,N-q 
q - appx . 1 to 3 

r * 0,1 m 
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Normalised Autocorrelation Function 


A A 

R„ / R , r-0, 1,2, ... ,m 


Note that: -1 ^ R f / R Q ^ 1 


Spectral Density Function 


m-1 


G' - 2H1 (R„ + 2 2E. K Cos (nrk/m) +(-l) k R m ) 

k o r in 

k * 0,1 m 


Inverse Filter 




G k t (2q+l) |H|V 


• (2q+l)^ 


k ■ 1 , 2 , . . . ,m 
2 


G; - CJJ (2q+l) 


Aperture Compensation 


r k ti a/2Hm 

G k l Sin (En a72HmT J 


k - 1,2 m 

0 


Smooth Power Spectral Density Function 

0.5 (G q + G l ) 

°' 5 < G m-l + G m> 

>k - °- 5 G k + °- 25 < G k-l + G k+1> 

k * 1 , 2, . . . ,m-l 


A 

G. 


A 

G 


m 


A 

G, 
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Typical Values 


5mm x 20mm . 

(4x) x (10x) - 4 Ox 

12.5 M x 500 m . 

«... nwcuai QUC DlZe 

System Magnification 

50mm .... 

.... Luecave Aperture 

lOrnm'* .... 

• • • • "‘ax . bpatial Frequency 

57. 

.... to (band width) 

2 Omm 

• • • • s * 0 «. rot* 

X - 2 x 10 3 

• • • • required Data Length (X) 

Zx 

• • . . « no. or data points 

104A 


1200 

• • . . OX “1 * 

^ f max ([2) (50mm* 1 

nnr 

10 

• • • • uur 

m - 1 1 

Sx • A w . Olmm 


) 

c TQmm"* 
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TYPE SPLIST 

FILE SPLIST 

SP1A-RLDR/H FOR SPi 

SP2A-RLDR/A FOR SP2 

SPLA-RLDR/n FOR SPL 

SPHA-RLDR/H FOR SPH 

SP23A-RLDR/f1 FOR SP23 

SP24A-RLDR/H FOR SP24 

SP1-HAIN WITHOUT FILTERS 

SP2-HAIN WITH FILTERS 

SPL-HAIN UITH LOU FREQ. FILTER 

SP3-PTDT 

SP4-SKPFL 

SPS-RDFL 

SP6-HDR 

SP7-C0RD 

SP8-RFP1T 

SPS'GBYTE 

SP10-HEANS 

SP11-SIGWA 

SP12-AUT0C0R 

SP13" (DELETED ) 

SP14-HIST 

SP15-DCLS 
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SP17-NRKLZ 

SP18-ZHEAN 

SP19-PSDF 

SP20-SPDF 

SP21-INF 

SP22-AC 

SP23-LFF 

SP24-HFF 
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ro 


SSS/TSl SP3 SPA SP5 SP| |P7 SP. SP. SPl. SPU 
SP16 SP17 SP18 SP19 SPSS SPS® FORT. LB 


SP12 SP14 SP1S * 


TYPE SP2A 
Rtr?/n SPS sp3 
SP16 SP17 SP 1 8 
R 

TYPE SPLA 
RLDR/fl SPL SP3 
SP16 SP17 SP18 

R 

TYPE SPHA 
RLDR/H SPH SP3 
SPl 6 SP17 SP18 


SPA SPS SP6 SP7 SP8 
SP19 SPSi SP22 SPS® 

SP4 SPS SP6 SP7 SP8 
SP19 SP31 SPS5 SPS® 


SP4 SPS SP6 SP7 SP8 
SP 19 SP21 SPSS SPS® 


SP9 SP10 
FORT. LB 


SP9 SP1® 
FORT. LB 


SP9 SPl® 
FORT. LB 


R 

TYPE SP23A 

RLL‘R/M SP23 FORT. LB 

R 

TYPE SP24A 

RLPR/H SP24 FORT. LB 


R 


SP11 SP12 SPl 4 SP1S ~ 

SPll SPIS SPl 4 SP15 

SPl 1 SPl 2 SP14 SP15 * 


A 


ooo o o oooooooooooooo 


TYPE SP1 

FILE SP1-MAIN WITHOUT FILTERS 
PROGRAMMER NAME: S.P. SINGH 
REVISION : 2 
DATEI4/19/79 

PURPOSE J THIS MODULE ACCEPTS USER INPUTS AND CALL ALL 
SUBSEQUENT MODULES 

RRH IS AUTOCORRELATION VALUES OF DATA 
GKC IS POWER SPECTRAL DENSITY FUNCTION 
GKH IS SMOOTHED POWER SPECTRAL DENSITY 
GKCP IS THE INVERSE FILTER FUNCTION 
GKC3P IS THE APERTURE COMPENSATION FUNCTION 
ITEMP, XTEMP ARE BUFFERS 
NU IS THE NUMBER OF WORDS PER RECORD 
NP IS THE NUMBER OF POINTS PER RECORD 
INTEGER P,Q 

DIMENSION RRHtlOl ),ITEMP< 1000 ), 

1 GKC<101 ),GKH(1©1) 

DIMENSION GKCP ( 101 ),GKC3P(101 ) 

DIMENSION NU ( 50 ) , NP ( 50 ) 

DIMENSION XTEMP( 1000 ) 

ACCEPT 'FILE NUMBER, NF -', NF 
ACCEPT 'DIFFUSE DENSITY FACTOR, DDF-' .DDF 
ACCEPT 'MAXIMUM LAG FACTOR, M-',M 
ACCEPT 'HEIGHT OF SCANNING SLIT, XL-', XL 
ACCEPT 'WIDTH OF SCANNING SLIT,XA-',XA 
ACCEPT 'LOU FREQUENCY FILTER PARAMETER, P- ', P 
ACCEPT 'HIGH FREQUENCY FILTER PARAMETER, Q- ', Q 
P-0 
Q-0 

CALL PTDT (ITEMP, NF,H. NR. NW.NP) 

CALL F0PEN(4, 'PAR' ) 

WRITE BINARY (4) NP.P.Q.NR 
CALL FCLOS ( 4 ) 

A-0 

B-3200 
NC- 100 

CALL HIST ( ITEHP, A, 3, NC , NR , NP ) 


CALL NEANS(ITEnP,XTENP,NP,NR,DDF ) 

C CALL SUAP< •SP23.SV , ,IER)j CALL LFF(XTENP,TENP2,TENP3,NP,P,NR ) 

C CALL SUAP< •SP24.SU*, IER)j CALL HFF (XTENP,TENP2,TEHP3,NP,Q # NR ) 

CALL SIGf1A(XTENP,NR,NP > 

CALL AUTOCOR ( XTERP, NP , NR, R, RRH, P, 0 ) 

CALL PSDF(RRH,N,H,XL,GICC> 

C CALL INF (GKC,GKCP,Q,R ) 

CALL AC1 ( GICC , GICC3P, H, H, XA ) 

CALL SPDF ( N, GKC3P, GKH ) 

CALL DELETE ( 'RD* ) 

CALL DELETE ( # NRD* ) 

CALL DELETE ( •ZRD* ) 

CALL DELETE ( 'ZP* ) 

CALL DELETE v 'TD* ) 

CALL DELETE ( 'PAR • ) 

STOP 

END 

R 


oooooooooooooo 


TYPE SP2 

FILE SP2-NAIN UITH FILTERS 
PROGRANNER NANE l S.P. SINGH 
REVISION'S 
DATE' 4/23/79 

PURPOSE* THIS NODULE ACCEPTS USER INPUTS AND CALL ALL 
SUBSEQUENT NODULES 

RRH IS AUTOCORRELATION UALUES OF DATA 
GKC IS POWER SPECTRAL DENSITY FUNCTION 
GKH IS SHOOTHED POWER SPECTRAL DENSITY 
GKCP IS THE INUERSE FILTER FUNCTION 
GKC3P IS THE APERTURE CONPENSATION FUNCTION 
ITENP, XTENP ARE BUFFERS 
NW IS THE NUNBER OF WORDS PER RECORD 
NP IS THE NUNBER OF POINTS PER RECORD 
INTEGER P,Q 

DINENSION RRH(101),ITENP<10C0>, 

1 GKC( 101 ),GKH( 101 ) 

DINENSION GKCP(101 >,GKC3P(101 ) 

DINENSION NW( 50 ), NP (50 ) 

DINENSION XTENP ( 1000 ) 

ACCEPT •FILE NUNBER, NF-',NF 
ACCEPT 'DIFFUSE DENSITY FACTOR, DDF- •, DDF 
ACCEPT • NAXINUN LAG FACTOR, N-*,N 
ACCEPT • HEIGHT OF SCANNING SLIT, XL-*, XL 
ACCEPT 'WIDTH OF SCANNING SLIT,XA-',XA 
ACCEPT 'LOW FREQUENCY FILTER PARANETER,P- ' ,P 
ACCEPT 'HIGH FREQUENCY FILTER PARANETER, Q- ' , Q 
CALL PTDT (ITENP, NF,H, NR, NU,NP) 

CALL FOPEN( 4, 'PAR' ) 

WRITE BINARY (4) NP,P,Q,NR 
CALL FCL0S( 4 ) 

A-0 

8-3200 
NC- 100 

CALL HIST ( ITENP , A , B , NC , NR, NP ; 

CALL NEANSdTENP, XTENP, NP, NR, DDF ) 

CALL SWAP( 'SP23.SU' , IER ); CALL LFF ( XTENP, TENP2.TENP3, NP,P, NR ) 


CALL SWAP ( * SP24 .SU* # IER); CALL MFF ! XTEHP , TENP2, TENP3 # NP, Q,NR ) 
CALL DELETE CZRD*) 

CALL RENAN ( •ZPP* , , ZRD*,IER) 

CALL SIGNA(XTENP,NR,NP> 

CALL AUTOCOR ( XTENP , NP, NR, N, RRH, P , 0 ) 

CAlL PSDF ( RRH,N, M, XL, GKC ) 

CALL INF !GKC,GKCP,Q,N) 

CALL AC < GKCP, GKC3P , N, H, XA ) 

CALL SPDF ! N, GKC3P, GKH ) 

CALL DELETE< *RD* ) 

CALL DELETE ( ‘NRD* ) 

CALL DELETE ( * ZRD* ) 

CALL DELETE! 'ZP* ) 

CALL DELETE! *TD* ) 

CALL DELETE! 'PAR* ) 

STOP 

END 


R 


oooooooooooooo 


TVPE SPL 

FILE SPL-MAIN UITH LOU FREQUENCY FILTER 
PROGRAMMER NAME: S.P. SIHGH 
REUISION : © 

DATE M/19/79 

PURPOSE: this module accepts user inputs AND call all 

SUBSEQUENT MODULES 

RRH IS AUTOCORRELATION UALUES OF DATA 
GKC IS POUER SPECTRAL DENSITY FUNCTION 
GKH IS SMOOTHED POUER SPECTRAL DENSITY 
GKCP IS THE INUERSE FILTER FUNCTION 
GKC3P IS THE APERTURE COMPENSATION FUNCTION 
T TEMP,XTEMP ARE BUFFERS 
.*U IS THE NUMBER OF UORDS PER RECORD 
NP IS THE NUMBER OF POINTS PER RECORD 
INTEGER P,Q 

DIMENSION RRH(101).ITEMP(1000), 

1 GKC( 101 ),GKH( 101 ) 

DIMENSION GKCP( 101 ),GKC3P( 101 ) 

DIMENSION NU ( 50 ) , NP i 50 J 
DIMENSION XTEMP(1000) 

ACCEPT 'FILE NUMBER, NF-',NF 
ACCEPT * DIFFUSE DENSITY FACTOR, DDF -', DDF 
ACCEPT 'MAXIMUM LAG FACTOR, M-*,M 
ACCEPT 'HEIGHT OF SCANNING SLIT. XL**, XL 
ACCEPT 'UIDTH OF SCANNING SLIT,XA-',XA 
ACCEPT 'LOU FREQUENCY FILTER PARAMETER, P- * »P 
C ACCEPT 'HIGH FREQUENCY FILTER PARAMETER, Q- • ,Q 

Q-0 

CALL PTDT (ITEMP,NF,H,NR,NU,NP) 

CALL FOPEN C 4 , 'PAR* ) 

URITE BINARY (4) NP , P , 0 , NR 
CALL FCLOS ( 4 ) 

A-0 

8-3200 

NC-100 

CALL HIST ( I TEMP, A, B,NC,NR,NP ) 

CALL MEANS (ITEMP,XTEMP.NP,NP f DDF ) 


CALL SUAP( •SP23.SU* , I ER ); CALL LFF C XTENP , TENP2, TENP1, NP, P, NR ) 
CALL DELETE ( 'ZRD* ) 

CALL RENAN ( 'ZP* , *ZRD* , IER > 

CALL SUAP ( * SP24 . SO * , IER ) ; CALL HFF ( XTENP, TENP2, TENP3, NP, 0 . NR ) 
CALL SIGN A C XTENP, NR , NP ) 

CALL AUTOCOR ( XTENP, NP, NR, N, RRH, P, Q ) 

CALL PSDF ( RRH, N,H,XL,GICC) 

CALL INF(GKC,GKCP,Q,N) 

CALL AC1 ( GKC,GKC3P, N,H, XA ) 

CALL SPDF(N,GKC3P,GKH> 

CALL DELETE ( 'RD* ) 

CALL DELETE! 'NRD* ) 

CALL DELETE ( •"RD" ) 

CALL DELETE ( *Z P* ) 

CALL DELETE! 'TD') 

CALL DELETE ( •PAR* ) 

STOP 

END 


oooooooooooooo 


TYPE SPH 

FILE SPH-MAIN WITH HIGH FREQ. FILTER ONLY 
PROGRAMMER NAME* S.P. SINGH 
REUISIONI0 
DATE : 4/24/79 

PURPOSE: THIS MODULE ACCEPTS USER INPUTS AND CALL ALL 
SUBSEQUENT MODULES 

RRH IS AUTOCORRELATION UALUES OF DATA 
GKC IS POUER SPECTRAL DENSITY FUNCTION 
GKH IS SMOOTHED POUER SPECTRAL DENSITY 
GICCP IS THE INUERSE FILTER FUNCTION 
GKC3P IS THE APERTURE COMPENSATION FUNCTION 
ITEMP, XTEMP ARE BUFFERS 
NU IS THE NUMBER OF UORDS PER RECORD 
NP IS THE NUMBER OF POINTS PER RECORD 
INTEGER P, Q 

DIMENSION RRH( 101 ), I TEMPI 1000 ) , 

1 GKC ( 101 ), GKH( 101 ) 

DIMENSION GKCP (101 ) , GKC3P ( 101 ) 

DIMENSION NU ( 50 ) , NP ( 50 > 

DIMENSION XTEMP< 1000 ) 

ACCEPT 'FILE NUMBER, NF- ■ ,NF 
ACCEPT 'DIFFUSE DENSITY FACTOR , DDF •', DDF 
ACCEPT 'MAXIMUM LAG FACTOR, M-',M 
ACCEPT 'HEIGHT OF SCANNING SLIT, XL-'. XL 
ACCEPT 'UIDTH OF SCANNING SLIT.XA-'.XA 
C ACCEPT 'LOU FREQUENCY FILTER PARAMETER, P- ' ,P 

P-0 

ACCEPT 'HIGH FREQUENCY FILTER PARAMETER, Q- ' ,Q 
CALL PTDT ( I7EMP,NF,H,NR,NU,NP ) 

CALL FOPEN( 4, 'PAR' ) 

URITE BINARY (4) NP,P,Q,NR 
CALL FCLOS( 4 ) 

A-0 

B-3200 
NC- 100 

CALL HIST ( ITEMP, A, B,NC, NR,NP ) 

CALL MEANS( ITEMP, XTEMP,NP, NR, DDF ) 


CALL SUAPCSP23. SUMER), CALL LFF(XTEHP,TENP2,TEHP3.NP,P,NR ) 
CALL RENAN ( *ZRD* , *ZP* , IER ) 

CALL SUAPCSP24.su*, IER), CALL HFF ( XTEHP, TENP2, TENP3, NP , Q , NR ) 
CALL DELETE ( *ZRD* ) 

CALL RENAfl ( * ZPP * , * ZRD * , IER ) 

CALL SI GNA( XTEHP, NR , NP ) 

CALL AUTOCOR ( XTENP, NP, NR, N, RRH, P, Q ) 

CALL PSDF ( RRH, fl, H, XL , GKC ) 

CALL INF (GKC,GKCP,Q,H) 

CALL AC(GKCP,GKC3P,n,H,XA) 

CALL SPDF ( N, GKC3P, GICH ) 

CALL DELETE( *RD* ) 

CALL DELETE( *NRD * ) 

CALL DELETE( *ZRD* ) 

CALL DELETE( *ZP* ) 

CALL DELETE ( *TD* ) 

CALL DELETE( *PAR* ) 

STOP 

END 


*VPE SP3 


C 

C 

C 

C 

C 

C 


ie 


o!i £ Sp 3-PTDT 

Revision^ NA " E,S - P * sing h 

DATE * 1 <'26/70 
d UP» 0S£ i TlJTC ma k 

S5S.ap«Lr 

c &wv£iy ' ,ER) 

8ft H> 

TVPE *Rau DATc P cr R ' fiUtNP) 

RETURN A STORED on DISK • 

END 


o o o o o o 


TYPE SP4 


10 


ae 


FILE SP4-SKPFL 

PROGRAMMER NANEiS.P. SINGH 
REVISION : 0 
DATE J 1/30/79 

SPEC?F!ED T ^ S T^ D UsiR SltIPS ™ E " U " 8EP 0F F “« 

COUP I LER NOSTACK 
SUBROUTINE SKPFL (INPBUF.NF) 

DIMENSION INPBUF ( 1000 ) 

KFL-1 

?? L l'lciI? I 2r 1 ?',®Ai N ?! UF ' ISTPT ’ IER - IC"T > 

lr iICNT.NE.l) GO TO 10 
KFL-KFL+1 

IF (KFL.EQ.NF) GO TO 20 

GO TO 10 

CONTINUE 

nfi-nf-i 

RE?5rA NU " BER 0,7 FILES SKIPP «".NF1 
END 


END 


TYPE SP5 


C 

C 

c 

c 


c 

c 

c 

c 


c 

c 

c 


le 

20 


C 

C 


100 


FILE SPS-RDFL 

PROGRAMMER NAMEJS.P. SINGH 

REVISION » 0 

DATE*2'2/79 

SUBROUTINE 0 RDFL*( INPBUF # NR,NU ) 

rd?l°rIads one file frgh tape into inpbuf 

AND WRITES IT IN DISK ?ILE TD * 

ICNT IS NO. OF WORDS IN A RECORD 
DIMENSION NW(S0).INP?UF(1000) 


NR IS NO. OF RECORDS 11M1M 

NU IS ARRAY CONTAINING NO. OF WORDS 
TYPE 'FILE IS BEING READ FROM TAPE 

NR ■ 0 

call fopen(2,*td*#2000> 

TYPE 'FILE TD OPENED* 

DO 10 I" 1 #50 
NUCI )-0 
CONTINUE 
DO 20 1-1,1000 
INPBUFCI )«0 

CALL 1 MTD 10(13,0# INPBUF, 1ST, IER. ICNT > 
TYPE "NO OF SAMPLES- *# INPBUF ( 2t ) 

TYPE *RDFL ICNT-*,ICNT 
IF C ICNT .EQ. 1 ) GO TO 100 
NR-NR+l 
NW(NR)-ICNT 

CALL WRITRC2.NR, INPBUF, l.IER) 

GO TO 10 

TYPE 1 'NUMBER OF RECORDS READ- ' * Stttpn 
TYPE 'FILE READ FROM TAPE AND WRITTEN 

CALL FCL0S(2> 

RETURN 


IN EACH RECORD 


ON DISK* 



oooooo 


TYPE SP6 

FILE SP6-HDR 

PROGRAMMER NAMEtS.P. SIMGH 
REUISI0N*0 
DATE *2/7/79 

PURPOSE* THIS MODULE READS THE TAPE DATA FROM FILE 
«TD', PRINTS HEADER INFORMATION, AND COMPUTES H-DELTA X 
COMPILER NOSTACK 
SUBROUTINE HDR( I TEMP, IREC,H ) 

DIMENSION ITEMP(1000) 

DOUBLE PRECISION IUD, IUD1 , IUD2, I I 
CALL F0PEN(2, 'TD* ,2000 ) 

CALL READR(2, IREC, I TEMP, 1 , IER ) 

WRITE (10,10) ( I TEMP (J ),J a 2,21 ) 

10 FORMAT (1X,20A2) 

TYPE 'NUMBER OF SAMPLES- ', ITEMP( 22 ) 

CALL CORD( ITEMP,23, IUD ) 

IX- IUD 

TYPE 'X COORDINATE (MICRONS )•', IX 
CALL CORD ( ITEMP, 25, IUD ) 

I Y- IUD 

TYPE 'Y COORDINATE (MICRONS IY 
CALL CORD( ITEMP, 27, IUD) 

IX- IUD 

TYPE 'DELTA X (H) (MICRONS )•', IX 
H- IX 

IF (H.LT.0) H--H 
CALL FCL0S(2 ) 

RETURN 
END 


O O O O O O J o o 


TVPE SP7 

FILE SP7-C0RD 

PROGRAMMER NAME: S.P. SINGH 
REVISION: 0 
DATE: 2/14/79 

PURPOSE: THIS NODULE DECODES X,Y COORDINATES AND 
DELTA X FRON 2 INPUT UORDS 
COMPILER NOSTACIC 
SUBROUTINE CORD< ITEMP, J, IUD ) 

DIMENSION ITENPC1000) 

DOUBLE PRECISION IUD, IUD1 , IUD2, 1 1 
DOUBLE PRECISION I J, IK 
TYPE 'CORD UAS CALLED* 

K-1.0D0 
IJ-32768.0D0 
11-65536. 0D0 
IXl-ITEMP(J) 

IX2-ITENPC J+1 ) 

TYPE *1X1-*, 1X1 
TYPE * 1X2- * , 1X2 
IF (IX1.LT.0) GO TO 10 
IF CIX2.LT. 0) GO TO 20 
IUD1-IX1 
IUD2-IX2 

IUD- IUD1XI I+IUD2 
GO TO 50 
20 CALL ICLRC 1X2.15) 

C TYPE * 1X2 AFTER ICLR*,IX2 

IUD1-IX1 
IUD2-IX2 
IUD2- IUD2+I J 
IUD - IUDHI I + IUD2 
GO TO 50 

10 IX3--IX1-1 

1X4 — 1X2-1 

IFCIX4.LT. 0) GO TO 30 

IUD1-IX3 

IUD2-IX4 


IUD-IUD1SII+IUD2 
IUD'IUD+IK 
IUD--IUD 
GO TO 50 

CALL ICLR(IX4,15) 
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TYPE SP8 

FILE SP8-RFHT 

PROGRAMMER NAME* S.P. SINGH 
REVISION: 0 
DATE: 2/12/79 

PURPOSE: THIS MODULE REFORMATS EACH RECORD TO THE 
16 BIT INTEGER FORMAT OF THE NOUA 1200 AND 
PLACES IT IN FILE *RD* 

COMPILER NOSTACIC 

SUBROUTINE RENT ( ITEMP, NR, NU, NP ) 

DIMENSION ITENP(1000>,XT<300) 

DIMENSION NP(50),NU(5O> 

COMMON JTEMP(1000) 

EQUIUALENCE < ITEMP, JTEHP > 

CALL F0PEN(2 # *TD* ,2000 ) 

CALL FOPEN <3, *RD* ,2000 ) 

TYPE * DATA IS BEING FORMATTED* 

DO 20 JR- 1 , NR 
DO 10 1-1,1000 
JTEMPCI )-0 
C ITEMP(I)-0 

1* CONTINUE 

CALL READR<2. JR, ITEMP, 1, IER) 

NUl-HU(JP) 

C TYPE * NU1 • ' , NU1 

DO 30 J-29.NU1 

C TYPE ' ABOUT TO CALL GBYTt * 

INP-ITEMP(J) 

CALL GBYTE CINP,IOU,0,12) 

C TYPE *GBVTE UAS CALLED* 

JTEMP( J-28 )*IOU 
3e CONTINUE 

C TYPE 'ABOUT TC CALL URITR* 

CALL URITR ( 3, JR, JTEHP, 1 , IER ) 

C TYPE ‘URITR CALLED* 

NP(JR)-NU< JRJ-28 
20 CONTINUE 

CALL READR C 3. 1 , JTEHP, 1 , IER ) 


II- 




49 

59 


R 


HPl-NP(l) 

IF iNPl.GT.399) NPl-399 
DO 46 I • 1 # NPl 
JK* JTEWP( I )/8 
XT( I )-JK/100. 

CONTINUE 


URITE (19,59) (XT ( I ), I 
FORMAT ( 1 9 ( 1X.F5.2 ) ) 
TVPE * DATA FORMATTING 
CALL FCLOS ( 2 ) 

CALL FCLOS( 3 ) 

RETURN 

END 


■1.NP1 ) 
COMPLETED* 


oooooooooo 


I 

rc 
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TYPE SP9 

FILE SP9-GBYTE 
PROGRAMMER NAME * S.P. SINGH 
REUISION » 0 
DATE* 2/16/79 

PURPOSE * THIS MODULE EXTRACTS BIT PATTERNS FROM 
A 16 BIT INTEGER UORD. 

IPK IS THE 16 BIT UORD 

IUPK IS THE BIT PATTERN REMOUED FROM IPK 
NOFF IS THE STARTING BIT POSITION IN IPK 
NBTS IS THE NUMBER OF BITS TO BE REMOUED FROM IPK 
COMPILER NOSTACK 

SUBROUTINE GBYTE( IPK, IUPK, NOFF, NBTS > 

INTEGER MASK( lS).LOC(lS) 

DATA LOC/1K,2K,4K,10K,20K,4OK,100K,200K,400K. 
t 1000K , 2000K , 4000K , 10000K , 20000K , 40000K/ 

DATA MASK/77777K , 37777K , 1 7777K , 7777K . 3777K , 1 777K , 
* 777K , 377K , 177K, 77K, 37K , 17K, 7K, 3K, IK/ 

IUPK • IPK . AND . MASK ( 1 ) 

IF (NOFF . LT . 2 ) GO TO 5 
IUPK* IUPK . AND . MASK ( NOFF ) 

5 NSHFT • 16- ( NOFF +NBTS ) 

IF(NSHFT .EQ.0.OR.NSHFT.EQ. 15) GO TO 12 
IUPK- IUPK/ (2**NSHFT) 

12 CONTINUE 

IF(NSHFT.EQ.IS) IUPK-0 

IF ( NOFF . GT . 0 .OR . IPK . GE . 0 ) RETURN 

IUPK* IUPK. OR. LOC (NBTS) 

RETURN 

END 

R 


o o o o o o - 


TVPE 


9 


SP10 

FILE SP10-NEANS 
PROGRANNER HARE * S.P. SINGH 
REUISIONt 1 


DATE * 3/19/79 

PURPOSE* THIS NODULE CALLS OTHER NODULES TO CONPUTE NEAN OF 
COMPILES SSsTS2r ALI2E T ° DIFFUSE UNITS AND ZERO NEAN 


SUBROUT I NE NE ANS ( ITENP , XTENP , NP . NR , DDF ) 
SS10 IS NEANS 
DINENSION ITENP(ieee) 

DINENSION XTENP(ieee) 

DINENSION NPC 50 ) 

ICHR-3 

ICHN-4 

ICHZ-5 


CALL AUGC ITENP, ICHR, NR, NP,RAUN> 

CALL NRRLZ< ITENP, XTENP, DDF , ICHR, NR , NP, ICHN, RAUN ) 
DDN-DDF 

TYPE *NEAN OF NORNALIZED DATA- • , DDN 
CALL ZNEAN( XTENP, ICHN. ICHZ, NR , NP, DDN ) 

RETURN 

END 


1 


type SP11 


I 

to 

to 


c 

c 

c 

c 

c 

c 


c 

c 


se 

le 

c 

c 


FILE SPll-SIGMA 
PROGRAMMER NAME * S.P. SINGH 
REVISION* 0 
DATE* 3/13/79 

mi ■ %"r bean souar£ uftLuE 

COMPILER NOSTACJC 
SUBROUTINE SIGMA (TEMP, NR, HP) 

ICHZ N | I0N TEHP( 1600 ),NP(50 ) 

CALL FOPEN( ICHZ, *ZRD*,‘<0OC ) 

SSQ*0 . 0 
NPI-0 

DO 10 I ■ 1 , NR 

NPi^NPa? R(ICHZ,I ' TEnP#1 ' IER) 


MEAN SQUARE 


THIS COMPUTES THE 
DO 20 J-1,NPI 

SSQ-SSQ+(TEMP( J )»t2 ) 

CONTINUE 

INP-lNP-fNPI 

CONTINUE 

XSQ-SSQ/INP 

cTi S cS2? PUTES THE R00T HE:an square 
SIG-SQRT(XSQ«INP/(INP-1 )) 

type 'MEAN SQUARE UALUE • • XSG 
TVPE 'ROOT MEAN SQUARE UALUE • ' SIG 
CALL FCLOS(ICHZ) # 1G 

TYPE 'SIGMA UAS CALLED* 

return 

END 


ooooooooo 


TYPE SP12 

FILE SP12-AUT0C0R 
PROGRAHHER NAHE* S.P. SINGH 
REVISION* l 
DATE* 4/4/79 

PURPOSE* THIS NODULE CONPUTES THE AUTOCORRELATION AND 
N0RNALI2ED AUTOCORRELATION FUNCTION OF ZNEAN DATA. 

RAT IS THE AUTOCOR VALUES 
RATN IS THE NORNALI2ED AUTOCOR VAULES 
TEHP, XT ARE BUFFERS 
CONPILER NOSTACK 

SUBROUT I NE AUTOCOR ( TENP , HP , NR , N , RAT , P , Q ) 

INTEGER P,Q 
DIHENSION NP(50 ) 

DIHENSION RATNt 101 ),RAT < 101 ) 

DIHENSION XT( 101 ) 

DIHENSION TEHP( 1000 ) 

TYPE 'AUTOCORRELATION IS BEING COHPUTED' 

ICHZ-S 

CALL FOPEN(ICKZ,'ZRD',4000) 

DO 5 I ■ 1 # 1000 
5 TEHP( I )»0 

DO 10 1-1,101 
RAT (I )-0 . 0 
XT( I ) *0 

10 CONTINUE 

HR-H41 
NR1-NRT1 
INP-0 

DO 20 I ■ 1 # NR1 

IF( I .EQ.NR1 ) GO TO 90 

CALL RE ADR ( ICHZ, I, TEHP, 1 , IER ) 

INP- INP+NP ( I ) 

NP1-NPCI ) 

DO 30 13-1, NP1 

TEHP( NP1 ♦H-1 3+1 >-TEHP(NPl-I3*l ) 

30 CONTINUE 

90 CONTINUE 


11-24 


40 

50 


70 

60 

80 

20 

85 

100 


1000 

2000 

R 


DO 40 14-1, H 

TEf1P(l4)-XTU4 ) 

CONTINUE 

DO 50 15-1 , n 

XT( 15 )-TENP(NPl + I5 ) 

CONTINUE 
DO 60 J-1,NR 
SUN-0. 0 
JR - J- 1 

DO 70 17-1, NP1 

CONTINUE TE ” P< 1 7 >tTENP ( I7 * JR ’ 


RAT( J )-RAT( J )-fSUN 
CONTINUE 

IF( I .NE.NR ) GO TO 20 
DO 80 18-1,1000 
TENP ( 18 ) -0 . 0 
CONTINUE 
CONTINUE 


RAT( J ) - RAT ( J )/ ( INP-J4 

continue 


-C2*(P*0>)) 


DO 100 J - 1 , NR 

RATN( J )-RAT( J )/RAT( 1 ) 

CONTINUE 

TVPE • R 

DO 2000 J-i.NR 

IR-J-1 


FORNAT ( lx! 1 3^ 8X ^1 ( 8X # F 1 1 N 4 
CONTINUE ' 1J ' 8 *' E11 - 4 ' S *'E11.4 


CALL FCLOS(ICHZ) 

RETURN 

END 


) 


NORN. 


t 

$ 

¥ 


R-HAT" 


r * 
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TVPE SP14 

FILE SP14-HIST 
PROGRAMMER NAME: S.P. SINGH 
REVISION* e 
DATE* 2/20/79 

PURPOSE* THIS MODULE COMPUTES A HISTOGRAM FOR THE DATA 
COMPILER NOSTACK 

SUBROUTINE HIST ( ITEMP,A, B,NC,NR,NP ) 

DIMENSION I TEMP ( 1000 ) , NP( 50 ) , NHST ( J 01 >,NST( 101 ) 

CALL FOPEN( 3, *RD* ,2000 ) 

TVPE ’HISTOGRAM BEING COMPUTED* 

Cl IS CLASS INTERVAL 
NC1-NC+1 
Cl • ( B-A )/NC 
DO 5 J-1.NC1 
MST( J )-0 
CONTINUE 

COMPUTE HISTOGRAM 
DO 10 J-l.NR 

CALL READR(3,J,ITEMP,1 # IER) 

NPl-NP(J) 

CALL DCLS(ITEMP,NP1,CI,NHST> 

DO 20 I ■ 1 , NCI 
MST(I )-MST(I ) + NHSTCI ) 

20 CONTINUE 

10 CONTINUE 

COMPUTE FIRST CLASS CONTAINING AT LEAST ONE POINT 

DO 30 >1.NC1 
IF (Mf r < I ).GE.l ) GO TO 40 
30 CONTINUE 

40 ININ * i 

C 

C COMPUTE LAST CLASS CONTAINING AT LEAST ONE POINT 
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TYPE SP15 

FILE SP1S-DCLS 
PROGRAMMER NAME* S.P. SINGH 
REUISION: e 
DATE i 2/26/79 

PURPOSE l THIS MODULE ASSIGNS A DATA UALUE TO A CLASS 
INTERVAL FOR COMPUTING HISTOGRAM. 

COMPILER NOSTACIC 

SUBROUTINE DCLS< ITEMP,NP1 .Cl .NHST ) 

DIMENSION I TEMP ( 1000 ) , NHST (101) 

DO 10 1-1,101 
NHST ( I )-0 
10 CONTINUE 

DO 20 I ■ 1 , NP 1 
J-IFIXC ITEHP( I )/CI ) 

NHST ( J+l >-NHST< J*1 ) + l 
20 CONTINUE 

RETURN 
END 
R 


o o o o o o o oo 


TYPE SPlfi 

FILE SP16-AUG 

PROGRAMMER NARE : S.P. SINGH 
REVISION * 8 
DATE* 3/6/79 

PURPOSE i THIS MODULE COMPUTES THE RAU MEAN OF THE DATA. 
COMPILER NOSTACK 

SUBROUTINE AUG( ITEhP, ICHR, NR, NP, RflUR ) 

SP16 IS AUG 
REAL ISUR 

DIMENSION ITERP(10e0).NP(50) 

CALL FOPEN( ICHR, 'RD* ,2880 ) 

INP-8 
ISUR-0 

SUH ALL THE POINTS IN THE FILE 
DO 10 I-l.NR 

CALL RE ADR < ICHR, I , I TEMP, 1 , IER ) 

NPI-NP(I ) 

DO 30 J» 1 , NFI 

isun-isun+iTERP(j) 

20 CONTINUE 

INP-INP+NPI 
10 CONTINUE 

RAUR-ISUR/INP 

TYPE * RAU MEAN - * , RAUR 

CALL FCLOS ( ICHR ) 

RETURN 

END 


R 
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TYPE SP17 

FILE SP17-NRHL2 
PROGRAHNER NAPE i S.P. SINGH 
REVISION: 1 
DATE: 3/19/79 

PURPOSE: THIS NODULE NORHALIZES THE FILE 
DATA BY HULTIPLYING THE DATA POINTS 
BY THE DIFFUSE DENSITY FACTOR AND CREATES A NEU FILE. 
CONPILER NOSTACK 

SUBROUTINE NPNLZ( ITEHP, XTENP, DDF , ICHR, NR. NP, ICHN.RAUN > 
C SP17 IS NRNLZ 

DIMENSION ITEHP< 1000 ),NP<50 > 

DINENSION XTENPC1000) 

CALL FOPEN( ICHN, • NRD* , 4000 1 
CALL FOPEN( ICHR , *RD ' , 2000 ) 

DO 10 I • 1 # NR 

CALL RE ADR ( ICHR , I , ITENP , 1 , IER ) 

NPI «-NP ( I ) 

DO 20 J • 1 , NP I 

XTENP ( J )- ITENP( J UDDF/RAUH 
20 CONTINUE 

CALL URITR( ICHN, I, XTENP, 1, IER ) 

10 CONTINUE 

CALL FCLOS(ICHR) 

CALL FCLOS(ICHN) 

RETURN 

END 

R 


I 
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TYPE SP18 


C 

C 

c 

c 

c 

c 

c 


S.P. SINGH 


FILE SPJ8-ZMEAN 
PROGRAMMER NAME : 

REUIS ION l 0 
DATE: 3/9/79 

PURPOSE: THIS MODULE NORMALIZES DATA BY SUBTRftAT TNr 

FR0 " ™ 

COMPILER NOSTACK 

SS18°iS^2MEAN^^^* I TEMP, ICHN, ICHZ,NR,NP, DDH ) 

DIMENSION I TEMP ( 1000 ).NP( 50 ) 

REAL ITEMP 

CALL FOPEN( ICHZ, ’ZRD* , 4000 ) 

CALL FOPENC ICHN, , NRD* , 4000 ) 

DO 10 I • 1 , NR 

CALLJEJDRC ICHN. I , ITEMP. 1. IER ) 


20 

10 


DO 20 J-l.NPI 

I TEMP (J )-ITEMP( J )-DDM 

CONTINUE 

CALL URITR(ICHZ,I,1TEMP.1 
CONTINUE 

CALL FCLOS(ICHN) 

CALL FCLOS(ICHZ) 

RETURN 

END 


IER) 


R 


ooo ooooo 


TVPE SP19 

FILE SP19-PSD r 
PROGRAMMER NAME x S.P. SINGH 
REVISION » e 
DATEx 3/19/79 

PURPOSE x THIS NODULE CONPUTES THE POWER SPECTRAL DENSITY FUNCTION. 
CONPILER NOSTACK 
SUBROUTINE PSDF ( RAT, H, H, XL. GKC > 

DIMENSION RAT ( 101 ).GKC( 101 ).GKH( 101 ) 

HR-HM 

COMPUTE THE POWER SPECTRAL DENSITY FUNCTION 

DO 10 KR-l.NR 
SUN-0.0 
NR2-NR-2 
DO 20 JR-1.NR2 

SUH-SUN+RAT ( JR-M )*COS< 3 . 1 4159* JR* ( KR-1 >/<HR-l )) 

20 CONTINUE 

GKC(KR)-2.*H*XL*(RAT< 1 )+2 . *SUH*RAT ( HR )*(-l )**(KR-1 ) ) 

10 CONTINUE 

RETURN 
END 
R 


type spae 


ro 


c 

c 

c 

c 

c 


c 

c 

c 


te 


10O 

2000 


file SP20-SPDF 
PROGRAMMER MAME: S.P. SINGH 
REVISION: 1 
DATE: 4/20/79 

coSp?l!r N0S?«cr DULE C0 " PL ' TES ™ E S " C0THEB oensitv function 

SUBROUTINE SPDF <M # GKC, GICH ) 

DIMENSION GKCdei ),G»:H(101 ) 

HR-M+i 

GKH(1 )-0.5*(G*CC(t )+GKC(2» 

GKH ( MR )«0 . St ( GKC ( HR- 1 )+GKC(MR)) 

MR1 -MR-1 

COMPUTE SMOOTHED DENSITY FUNCTION 
DO 10 ICR «2, MR 1 

CONTINUE® ' SSt< ° KC ‘ * R ~ 1 ’* G,a: < **♦ 1 > >♦« • S*«C < ICR ) 

™ P | * * GICC3POC ) 

DO 2000 J ■ 1 , MR 
IR-J-1 

™I!i ie,10e) IR »^CCJ).GICHCJ) 

RETURN 
END 


GH(K ) 


R 


oooooo 


TYPE SP21 

FILE 


C 

C 

C 


ie 


ie& 

2000 


SP21 • IHF 

PROGRAMMER NAME* S.P. SjNCH 

REVISION* 1 

PURPOSE » THIS flODULE PERFORMS 
Oh THE POWER SPECTRAL DENSITY 
COMPILER NOSTACK 
SUBROUTINE INF (GKC '5££o',?i? ! 
DIMENSION GICCC 101 ),GKCP( 1®1 

INTEGER Q 
PI -3. 14159 

MR-M+1 

COMPUTE THE FILTERED PSDF 

DO 10 Kl-l.HR 
K-Kl-1 

IF(K.GT.0> GO TO 5 
giccpcki )-G«ccna ) 

GO TO 10 


THE INUERSE 
FUNCTION. 


FILTER OPERATION 


IF(SIN( (K*PI ) ) «LT. 1 1^0E-06 ENCOUNTERED IN IHF* 

•THE DENOMINATOR OF SINE ^UNCTIOH KlOU ' /SIH( „ OT ),(8»W) 

GKCPCK1 ) - GKCCICi )tC C2t0^1 )*SIN( (K* 1 c 


GKCP 
/(21M) ) )**2 
CONTINUE 

TYPE ' K GKC ( K ) 

DO 2000 J • 1 # HR 

WRITE! 10 » 100 1 IR,GKC( J ),GKCP( J ) 

FORMAT (1X,I3#8X # E11»^»8X,E11.^ 1 

CONTINUE 

RETURN 

END 


GlCCPdC)' 
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TVPE SP22 

SP22*AC 

PROGRAHNER NAHE: S-P. SINGH 
REUISION: 1 
DATE: J/19'79 

PURPOSE: THIS NODULE PEREORNS THE APERTURE CONPENSATION 
OPERATION ON THE INUERSE FILTERED PSDF . 

CONPILER NOSTACK 

SUBROUTINE AC < GKCP, GKC3P, N , H, XA ) 

DINENSION GKCP( 101 ),GKC3PC lOn 
PI-3. 14159 
HR-N+1 

DO 10 K 1 • 1 , HR 
IC-K1-1 

IF ( K . GT . 0 ) GO TO 5 
GKC3PMC1 )-GICCP(Kl ) 

GO TO 10 

5 IF(SIN( <K*PI*XA)/(2*H*N) ).LT.l .0E-06) TVPE 

* * T HE DENONINATOR OF SINE FUNCTON BELOW 1.0E-O6 ENCOUNTERED IN AC* 

GKC3P< K 1 )-GKCP(Kl )*( ( (KtPIIXA )/<2IHIN ) )/SIN( (KJPItXA >/(2*H*H> > )t*2 

10 CONTINUE 

TVPE * K GKCP(K) GKC3P< K ) * 

DO 2000 J«1,HR 
IR-J-1 

WRITE (10, 100) IR,GKCP( J ) ,GKC3P( J ) 

10‘ FORHAT ( 1X # I3,8X,E11 . 4 , 8X , E 1 1 .4) 

2000 CONTINUE 

RETURN 
END 
R 


TYPE SP23 

C FILE SP23-LFF 

C PROGRAMMER NAME * S.P. SINGH 

C REUISION* 2 

C DATE* 4/26/79 

C PURPOSE* THIS NODULE PERFORMS THE LOU FREQUENCY OPERATION 

C ON THE Z-NORHALIZED DATA. 

CONPILER NOSTACK 

C SUBROUTINE LFF < TEMP 1 , TEMP2, TENP3. NPl , P, NR > 

DIMENSION TEMPI ( 1666 >,TEHP2< 1666 >,TENP3< 1660 ) 

INTEGER QU0,P,Q,P2,NP1(S6>,REC1,REC2,REC3 
TYPE * DATA IS BEING FILTERED FOR LOU FREQUENCY- 
CALL FOPEN^-ZRD'^eeO) 

CALL FOPENO, *ZP',4600 ) 

CALL F0PEN(4, -PAR- ) 

READ BINARY (4) NP1,P,Q,NR 
CALL FCL0SC4) 

C 

C CALCULATE t OR RECORDS TO COMPUTE FIRST SUM AND 

C FIND REC. t OF FIRST POINT 

C 

NP-NPKl ) 

NREC«(2*P*1 )/NP 

IF(M0D<2*P>1,NP).GT.C ' NREC-NRECM 
QUO-P/NP 

IF(NOD(P,NP).GT.0> QUO-QOO+1 
IF (QUO.EQ.l) GO TO 3 
JQUO-QUO-1 
DO 2 I ■ 1 , JQUO 

CALL READR ( 2, I , TEMPI , 1 , IER > 

CALL URITR (3,1, TEMPI , 1 , IER ) 

2 CONTINUE 

3 CONTINUE 
C 

C CALCULATE FIRST SUM (POINT P*l) 

C 

Hl-0 

SUM-6.6 


o o o ooni/i 


20 

10 


40 

C 


DO 10 I • 1 , NREC 

CALL READR(2,I, TEMPI, l.IER) 

DO 20 J-l.NP 
SUM-SUM+TEMPHJ) 

N1 • N1 ♦ 1 

IF(N!.EQ.(2*P+1 ) ) GO TO 5 

CONTINUE 

CONTINUE 

CONTINUE 

FIND LOCATION OF FIRST FOINT 

CALL REPDR(2,QU0,TtriPl,l,IER> 

LOCI "P+A-tNPT ( QUO-1 ) ) 

TEMPI ( LOCI )- TEMPI (LOCI >-(3UN/FL0AT(2*P+l >) 
CALL URITR( 3, QUO, TEMPI , 1 , IER ) 

P2-P+2 

N-NP*NR 

NMP-N-P 

DO 30 I-P2,NMP 


CALCULATE REC t OF PTS. I,I-P-1,I+P 


REC1-REC1+1 


REC2-REC2 i-1 


REC1-I/NP 
IF ( MOD ( I , NP > . GT . 0 ) 

RE J2- ( I-P-l )/NP 
IF(MOD(I-P-1,NP).GT.0) 

REC3-U + P )/NP 
IF(MOD(I+P,NP).GT.0 ) REC3-REC3+1 
IF (I.NF.P2) GO TO 40 
IREC1-REC1 
IREC2-REC2 
IREC3-REC3 
CALL READR 
CALL READR 
CALL READR 
CONTINUE 


<c,REC1,TEMP1 # 1.IER> 
(2,REC2,TEMP2, 1 , IER ) 
(2,REC3,TEMP3, 1, IER ) 


OO 1*1 ooooooo 





► 


C CALCULATE LOC. OF PTS. I.I-P-1,I*P IN RESPECTIVE RECORDS 

LOCI ■ I- ( NPt ( REC1-1 ) ) 

L0C2* ( I-P-l )-(NP*< REC2-1 ) ) 

L0C3* ( I *P )-(NPt (REC3-1 )) 

READ RESPECTIVE RECORDS TO OBTAIN I,I-P-1,I*P 
IF (IREC1.EQ.REC1 ) GO TO 5e 
CALL READR( 2, REC1 , TEMPI , 1 , IER ) 

IREC1-REC1 
50 CONTINUE 

IF ( IREC2.EQ.REC2) GO TO 60 
CALL READR(2,REC2,TEMP2, 1, IER) 

IREC2-REC2 
60 CONTINUE 

IF CIREC3.E0.REC3) GO TO 70 
CALL RE ADR ( 2, REC3, TEMP3, 1 , IER ) 

IREC3-REC3 
70 CONTINUE 

SUM-SUM-TEMP2( L0C2 )+TEMP3( LOC3 ) 

XT-TEMPKLOC1 )-(SUM/FL0AT(2*P*l ) ) 

CALL READR ( 3, REC1 , TENP1 , 1 , IER ) 

TEMPI ( LOCI ) -XT 

CALL URITR(3,REC1,TEMP1 # 1, IER ) 

0 CONTINUE 

IF (REC1.EQ.NR) GO TO 90 
REC1-REC1+1 
DO 80 I ■ REC 1 , NR 
CALL READR (2,1, TEMPI , l , IER ) 

CALL URITR (3, I, TEMPI, 1, IER) 

80 CONTINUE 

90 CONTINUE 

OUO-P/NP 

IF (MOD(P,NP ) .GT.0 ) QUO-QUO + 1 
IP-P-( (QUO-1 ' tNP ) 

CALL READRC3.QUO, TEMPI, 1, IER) 

TYPE m Z' m 

TYPE (TEMPI (I),I»IP,NP) 


11-38 


c- ZERO OUT BEGINNING AND LAST P 

C 

ICNT-0 
NP2-NP1 ( 1 ) 

NR2-NR 
DO 1 I - 1 , NR 

CALL READR<3,I,TENP1,1,IER) 

CALL RE A DR ( 3,NR2, TENP2, 1 , IER ) 

DO 280 J-1,NP2 

TEMPI CJ )-0.0 

TEMP2CJ)-0.0 

ICNT-ICNT+1 

IF(ICNT.EO.P) GO TO 380 
200 CONTINUE 

CALL URITR(3,I,TENP1,1,IER) 
CALL URITR(3,NR2,TEWP2,i,IER) 
NR2-NR2-1 
1 CONTINUE 

380 CALL URITRC3, 1,TEf1Pl, 1, IER ) 

CALL URITR(3,NR2,TENP2,1,IER) 


CALL PCLOS ( 2 ' 
CALL FCLOS ( 3 ) 
CALL BACK 
END 



POINTS FROM FILE ZP 


' 


i 
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TVPE SP24 

FILE SP24-HFF 

PROGRANNER NANE * S.P. SINGH 
REUISIONi 2 
DATE * 4/20/79 

PURPOSE i THIS NODULE PERFORHS THE HIGH FREOUENCV OPERATION ON 
THE LOU FREQUENCY DATA. 

CONPILER NOSTACK 

SUBROUTINE HFF ( TENP l , TENP2 , TEHP3, NP1 , P , NR ) 

DINENSION TEN?1( 1000 >,TENP2( 1000 ), TEHP3 ( 1000 ) 

INTEGER QU0,P,Q,P2, NP1 < 50 ) , REC 1 , REC2, REC3 
TYPE • DATA IS BEING FILTERED FOR HIGH FREQUENCY* 

CALL F OPEN (2, *2P* ,4000 ) 

CALL FOPEN< 3, *ZPP* , 4000 ) 

CALL F0PEN( 4, 'PAR* ) 

READ BINARY (4) NP1,Q,P,NR 
CALL FCLOS ( 4 ) 

CALCULATE * OF RECORDS TO CONPUTE FIRST SUN AND 
FIND REC. t OF FIRST POINT 

NP-NPK 1 ) 

NREC-(2*P + 1 )/NP 

IF ( NOD( 2*P+1 , NP ) . GT . 0 ) NREC-NREC+1 
QUO-P/NP 

IF ( NOD( P , NP ) . GT . 0 ) QUO-QUO+1 

CALCULATE FIRST SUN (POINT PM) 

Nl-0 
SUN-0.0 

DO 10 I • 1 # NREC 
CALL RE A DR (2, 1.TENPl, 1 # IER ) 

DO 20 J ■ 1 , NP 
SUN-SUN + TENPK J) 

N1-NH 1 

IF (N1 .EQ. <2*PM ) ) GO TO S 
20 CONTINUE 


o o o oomn 


10 CONTINUE 

CONTINUE 


40 

C 

C 


c 

c 

c 


FIND LOCATION OF FIRST POINT 

CALL READR (2, QUO, TEMPI, 1, IER ) 
LOCI -P+1 -<NP*< QUO- 1 ) ) 

TEMPI < LOCI ) ■ SUN/’FLOAT ( 2*P+ 1 ) 

CALL URI TR (3, QUO, TEflPl , 1 , IER ) 

P2-P+2 

N-NPtNR 

NNP-N-P 

DO 30 I-P2,NNP 


CALCULATE REC t OF PTS. I,I-P-i.I+P 


REC1-REC1+1 


REC2-REC2+1 


REC1-I/NP 
IF ( MOD( I,NP).GT.0) 

REC2- ( I-P- l )/NP 
IF(NOD( I-P-1,NP).GT.0) 

REC3-C I+P )/NP 
I F < MOD ( I +P, NP ).GT .0) REC3-REC3+1 
IF CI.NE.P2) GO TO 40 
IREC1-REC1 
IREC2-REC2 
IREC3-REC3 
CALL READR 
CALL READR 
CALL READR 
CONTINUE 


<2,REC1,TENP1,1, IER) 
<2,REC2,TEMP2,1,IER) 
(2»REC3,TENP3, 1,IER) 


CALCULATE LOC. OF PTS. 

LOC 1 » I - ( NPl ( REC1 - 1 ) ) 

LOC2* ( I-P-l )-(NP*(REC2-l ) ) 
L0C3- ( I+P )-(NP*(REC3-l ) ) 


I. I-P-l, I+P 


IN RESPECTIUE RECORDS 


READ RESPECTIUE RECORDS TO OBTAIN I, I-P-l, I+P 


ooooooooo 


r 






50 

60 

70 

30 


IF (IREC1 .EQ.REC1 ) GO TO 50 
CALL READR (2,REC1 .TEMPI , 1, IER > 

IREC1-REC1 

CONTINUE 

IF ( IREC2.EQ.REC2 ) GO TO 50 
CALL READR ( 2,REC2,TEHP2, 1 , IER ) 

IREC2-REC2 

CONTINUE 

IF ( IREC3.E0.REC3 ) GO TO 70 
CALL READR <2,REC3,TEMP3, 1, IER ) 

IREC3-REC3 

CONTINUE 

SUM-SUM-TEMP2(LOC2 >+TEMP3(LOC3 ) 

XT-SUM/FLOAT(2*P-*l ) 

CALL READR(3#REC1,TEMP1, 1, IER ) 

TEMPI ( LOCI >-XT 

CALL URITRC 3, REC1 , TEM°1 # 1 , IER ) 

CONTINUE 

QUO-O/NP 

IF ( MODCG, NP ) . GT . 0 ) QUO»QUO+l 

IP-Q-( (OUO-1 )*NP UP 

CALL READR ( 3, QUO, TEMPI , 1 , IER ) 

TVPE "Z"" 

TYPE ( TEMPI (I) # I"IP,NP) 

ZERO OUT BEGINNING AND LAST P+Q POINTS FROM FILE ZPP 


ICNT-0 
NP2-NPK1 ) 

NR1-NR 
DO 1 I • 1 # NR 

CALL READR (3,1, TEMP 1,1 # IER) 

CALL READR ( 3, NR1 , TEMP2, 1 , IER ) 

DO 200 J-1.NP2 

TEMPI ( J ) *0 . 0 

TEMP2( J )-e.0 

ICNT ■ ICNT+1 

IF ( ICNT.EO. (P+Q ) ) GO TO 300 




i 


4 


<*■ 


I 



200 CONTINUE 

v CALL URITRC3, !, TEMPI, 1 # IER) 

CALL URITR( 3, NR1 , TENP2, 1 , IER ) 
NR1-NR1-1 
1 CONTINUE 

300 CALL URITRO, I .TEMPI, 1, IER ) 

CALL L'RITRC 3, NR1 , TEMPS, 1 , IER ) 

C 
C 

CALL FCL0S(2 ) 

CALL FCLOSO) 

CALL BACK 
END 



I 






. 



TYPE spas 
C SP25-AC1 

C PROGRAMMER NAME * S.P. SINGH 

C REVISIONS 

C DATE 14/20/79 

C PURPOSE* THIS MODULE PERFORMS THE APERTURE COMPENSATION 

C OPERATION ON THE INVERSE FILTERED PSDF FOR SP1 AND SPL. 

COMPILER NOSTACIC 

SUBROUTINE AC1 < GICCP . GICC3P, H, H. XA ) 

DIMENSION GKCP<101 >,GKC3P(101 ) 

PI-3.14159 

MR-M-H 

DO 10 K1 ■ 1 , MR 
K-Kl-1 

IFUC.GT.0) GO TO 5 
GKC3P(K1 ) -GKCP ( K 1 ) 

GO TO 10 

5 IF ( SIN( < K*PI tXA )/< 21HIM ) ) . LT . 1 . 0E-06 ) TYPE 

t 'THE DENOMINATOR OF SINE FUNCTON BELOU 1.0E-06 ENCOUNTERCD IN AC* 
GKC3P (K1 ) -GKCP (K 1 )* ( < < K*PI *XA )/ ( 2*HtH ) >/SIN< < K*PI *XA )/ (2tH*M > ) )**2 
10 CONTINUE 

TVPE * K GKC(K) GICC3P<K>" 

DO 2000 J- 1 , MR 
IR-J-1 

WRITE ( 10, 100 ) IR,GKCP( J ),GKC3P( J ) 

100 FORMAT (IX, I3,8X.E11.4,8X,E11 .4 ) 

2000 CONTINUE 

RETURN 
END 



DIFFUSE PFNSITV FaCThR.DPF' 1.27 
NAXlNon FACTOR. 11*28 
hIIGmT OF SCAtWINi SLIT . XL* IOOO 
UIDTH OF jr.rf<MIHG ;IIT.xh-S? 
data preparation, coriftNcer 

(4JNB1 R OF FILLS -HIPFIP. * 

Far IS BUNG RIAL *B,N Tap* 
NU*S1P OF RECORDS PlAp* 



ROO* HF AM SQUARE oAtut - 9.2622 981 -| 

AUTOCORRELATION is icing confuted 


ft 

R-mAT 


NORN. R-MAT 

« 

• 68721 

-3 

C 1000C 


1 

I.4QSX 

-3 

0.72071 


2 

0 ■ 2693E 

-3 

0 79191 


] 

• . 171 If 

-3 

0.24901 


4 

O.ll^OC 

-3 

0 . 1 "’021 


S 

0.8441C 

-4 

0 12281 


6 

• 906 1C 

-4 

0. 1 31 91 


7 

0. 12021 

-3 

0.I74S1 


8 

0. 1SS6C 

-3 

0 .226*1 


9 

0. 1870E 


0.27211 


II 

0. 1820C 

-3 

0.26*81 


1 1 

0.1S19E 

-3 

0.22111 


12 

0. 1382F 

-3 

0.201 11 


13 

0. 132IC 

-3 

0.19221 


14 

0. I163C 

-3 

0.16921 


IS 

0. 11821 

-3 

0. 17181 


16 

0.13*11 

-3 

0. 189)1 


17 

*. 13921 

-3 

0.20261 


18 

0. MOSt 

-3 

0.20441 


19 

*. 17961 

*3 

0.20311 


20 

0 . 1 39S1 

-3 

0.2030C 


21 

0. I40SC 

-3 

0.204 91 


22 

e. 13271 

*3 

0.1 9301 


23 

*. U-2F 

-3 

0.17061 


24 

* . 103SC 

-3 

0.1S07E 


2S 

«.9isie 

-4 

0 13321 


26 

0.7764C 

-4 

* . 1 1 3M 


27 

0.8O14C 

4 

0 . 1 1661 


28 

0.9*311 

-4 

0.13721 


K 

<*CMf » 


UTC3PIK > 


e 

0.3133E 

3 

0.31371 

3 

t 

*.69’4f 

2 

*.69981 

2 

2 

2. b«bt 

2 

0.36761 

2 

3 

0. 57*21 


O.S847r 

2 

4 

0.41731 


0. 4J1U 

2 

s 

0.62811 

2 

0.66*81 

2 

6 

0.58831 

2 

0.63311 

2 

7 

• .33’Wt 

2 

0 4 3881 

2 

8 

0.26I6C 

2 

0.29831 

2 

9 

0 . ISfeSt 

2 

0 18771 

2 

10 

0.18341 

2 

0 22SS1 

2 

11 

0. 1 «bOC 

2 

0.12771 

2 

12 

0. 16S0C 

2 

0 22281 

2 

13 

0. 14961 

2 

0.21 321 

2 

14 

0. 11931 

2 

0. 16041 

2 

IS 

0 64*91 

1 

0.10341 

2 

16 

0. 882 it 

1 

0.1S261 

2 

17 

».7y<C 

1 

0. 13661 

2 

18 

0.80741 

1 

0.163*1 

2 

19 

0. 3621C 

1 

0 S002E 

1 

20 

0. 12351 

1 

0 78S71 

l 

21 

8.21941 

1 

0 59001 

1 

22 

0.288 1C 

1 

0 86461 

1 

2J 

e. i~b6t 

1 

0 59601 

1 

24 

0 18821 

1 

0 72061 

1 

2S 

0.12691 

1 

0 SS721 

1 

26 

*. I81SC 

1 

0 92 711 

1 

27 

0.11 tOf 

1 

0.69171 

1 

28 

0. 18871 

1 

O.7-*09C 

1 

K 

C«C3P<* 

i 

CM. 

K 1 

0 

0. 31331 

3 

0.19161 

3 
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SPW 

rite NIJNBfR.NF -5 
Difim* Ptrosm r^cro«.Ptf *i.27 
NAMHUB LmC »A.TuB.--£6 
mEICNT v<f SCANNING -IIT.»L-I002 
UlPTM Of x -hMf*. :>LIT.»a-2S 
hJuM FSfO««N4.v FILTER aaRm/Ic T£r . n- 1 
DaTA >.cr*“tr*.f D 

NUfiBER CO f ! L£ SflRFfD* * 

i III I> MINC Sl«»* FRON TmFI 
NUMBER Of OtCOSff R£aD* 

FILE •'far TATI AMP -S! T ’fN 0*» Dl! 

ISA 2 I ’ L f-l J'a?vT*SUP 1*' 1 1 i 
NuNIfR c APPl £ i - 261 

x cooeriNMiE «nic»ohs»* 

V COORDINATE t«|CSoMS »• * 

DELTA > . M> . HU RON* >• IS 

Data l , ££ 1 n »7 FORMATTED 


1.7b 1.79 


1.69 ISO I ■ 8 j 1.82 
1.85 1.82 1.84 1.87 


1. 81 t.7? 1.79 1.83 1.82 t-82 

1.66 1-82 1.81 1.84 1.82 1.84 

1.86 1.89 1.83 I -85 1.82 I . 88 

1.82 1.85 1 82 1.82 1.81 1-82 


1.79 1.62 1.81 

1.86 :.68 1 . 3 " 


1.8 1 I.B4 

1.79 1.79 


1.84 |.S4 1.89 IBS I.FO 1 .St 


1.84 1.34 1.86 


1.79 I. 81 


I. 82 1.80 1.79 1.79 |.8S 1.85 

J. 79 1.78 1.80 1.82 1 .89 1.84 

1.76 1.74 I.7S 1.78 1.79 1.78 


1.77 1.79 


1.78 1.78 1.78 

’..80 1.80 1.81 


1.82 1.62 1.89 1.63 1.84 1.81 

1.82 1.87 1.83 1.81 1.20 1.31 

1.81 I.7S 1.73 1.78 1.79 1.7S 

1.88 1.82 1.82 1.24 1.83 1 SO 

1.87 1.85 1.84 1.24 |.88 1.88 


1.83 *.81 

1.77 1.77 


1.82 1.84 


1.77 1.84 


1.88 1.88 
1.84 1.83 




Data 

neon 

IS OC INC FILTERED FOB 
Si-AR* <>*U,E • 0. 




BOO’ 

N£ AN SOUARf UALUf * 

0 2204 354 

>1 



OilTOCOOBf LAT10N 15 If INC c 

twin d 




a 

6 -MAT 

NOON. R-MA t 



0 

C.4B65E -1 

0. ioo«( 




• 

0.4I73C -3 

0.8578E 




2 

4. OSOf -3 

• 6269E 




3 

9 . £01 94 -3 

0.0 sot 




0 

0.14I5C -3 

0.29*9* 




5 

O.llSTf 3 

0 2377f 




6 

0.IJ33C -3 

0.23404 





0.13054 -3 

0 AIK 




8 

0.1524E -3 

0.313.-* 




9 

0. 166 Of -3 

0. 302 Sf 




to 

0. 16571 -3 

0.34#7f 




II 

0. 1 S52f -3 

0 320 # 




12 

0.1024f 1 

0.2527* 


1.82 

1.79 

13 

0 1336E -3 

O^’OSE 


1.79 

1.81 

10 

0 1295E -3 

0.2661E 


1.60 

1.85 

IS 

0.1329E -3 

0.2’33f 


1.81 

1 85 

1C 

0.I376C -3 

0 2829t 


1.80 

1.81 

:7 

0.1359C -3 

0 2876f 


1.85 

1 .82 

18 

O.I387f 3 

0.28SU 


1-82 

1 .85 

19 

0 . 1 381f -3 

0 28 394 


1.81 

1.86 

20 

0.13S71 -J 

0.27t9f 


1 83 

1.84 

21 

0. 1 303* -3 

0.2677E 


I.Tg 

1.81 

22 

0.I22OT -3 

0 251 £4 


l.7» 

1.8# 

23 

0.11 384 -3 

0.23404 


1.79 

1.78 

20 

O.I04Jf 3 

0.2I35C 


1.76 

1.77 

25 

f.'AVU -0 

0.I964C 


1.77 

1 .80 

26 

0.8979* 4 

0 18461 


1.84 

1.83 

27 

0.895X 4 

0.1840C 


1.80 

1.77 

28 

0.929-# -0 

0 19114 



1.79 1-83 1.79 1.78 


1.78 l.?6 

1 .81 1.65 184 

..85 

183 

1.81 

1 84 

1.60 

Tf# 

DO* r jnmTOR Of S 

INI 

1.77 1.-8 

1.79 1.77 1.81 

1.79 

1.76 

1.79 

1 .77 

1.80 

’Ml 

Dtr*07SNMT0« Of 

1.84 1.85 

1.87 1.84 1.75 

1.79 

S .82 

1.62 

1.78 

l.7fc 

TMt 

jfN fiNM -* c* s :»« 

1.78 1.78 

180 I.SO 1.80 

1.78 

1.79 

1.81 

1.79 

1 80 

1 

c». <* < 


1.77 1.75 

1.79 1.81 1.77 

1.77 

1 .78 

1.74 

1.70 

1.77 

0 

0.31I3C 

3 

1.79 







1 

0.65I8E 

2 

DATA FOOnATTlNG CONftf TED 






2 

0. JO! 24 

2 

RAW data S-OB4D LIST. 






9 

0.072LE 

£ 

MISTOCBAN IE INC CCnF>trED 






4 

0 3*554 

2 








5 

0.42204 

2 

TOT#iL no. Of 

CLASSES* 100 






6 

0. 38294 

2 

CLASS NO. 

NO. Of DATA POINTS 





7 

0.2415c 

2 

4| 

2 






8 

0. 13685 

2 

42 

2 






9 

0.69N2C 

1 

43 

2 






10 

0.7I0-* 

1 

44 

34 






11 

0 .521 4f 

1 

4$ 

206 






12 

0.4250* 

1 

46 

3«9 






13 

0.264 It 

1 

47 

206 






10 

0.2779* 

1 

48 

41 






15 

0 12704 

1 

49 

1 






16 

0 . l * 47 f 

1 

RAU KAN • 

0.105498E 4 






17 

0.'»2 7 * 

0 

Nf an Of NCBNAL I Zf D DATA. 

.12-0004 1 




18 

0 831 If 

0 


T Mt Pf NONIfAil A R 16 S|N£ ( Uf* T IuN Bf 19« I .01 

** 54 -.rnINAT’fl Of Sl-*i F UNC ~'."f IE LOW 1 Of 

’►* PENCJIINaT B CO S|Nt » UNC' ION BELOW I.OC 

*HC I'f' SIlU - Of 1N£ i UNC ' I liN if LOW I Of 

•Mt Of xnihat:s o» siNi FiC4r':«44 klow i.of 

•me PE non: mat o of siNf function iclcoj i.ot 

•me wncninatoo co sir# Fi*«c’ico4 suow i »l 


i«CP«i * 
0.3I18E 3 

2 .65 ? Jf 2 
0.ill*f 2 
8.S0»7f 2 

0.396 It 2 
0-S236C 2 

O.S244f 2 

• 3729C 2 

O 24J8F 2 

O 14584 2 

• saiir 2 

O. 1 '•21 2 

O.lS'Sf 2 

». r#n 2 

0. 25«1E 2 

f.iw 2 

1. 30*5* 2 

0.59*71 2 

0. «277| 3 


06 ENCC JfTfBfD IN INF 
06 £fw.0O6'fO*D IN INf 
4Nt tNOOUNTIKD IN I*# 
Oe CNC0«NT|PE6 ;n JNF 
06 ENCCc-NTfOfD IN II# 
06 fNCGUN’.CBf t IN II# 

«6 fN. vXJNT£A£f ;n ># 
-06 06 OUNTEBEfr IN :*# 
06 ENcOuNTfBft IN ir# 
of l r«r Ponte of r In :»# 


II -4!) 


19 

9.46396 

• 

9.19156 

4 

29 

9.77956 

• 

9.11596 

3 

21 

9 SI19£ 

9 

9.26816 

2 

22 

9.74626 

9 

9.21146 

2 

23 

9.54566 

9 

9.19216 

2 

24 

9 . 7 4 W* 

9 

9. 19466 

2 

25 

9.66766 

9 

9.76246 

1 

26 

9 . 99296 

9 

9.89916 

1 

27 

9 . 60696 

9 

9 55946 

1 

28 

9.47396 

9 

9 42656 

1 

K 

* 


Cl C >■ 4 1 


9 

9. 3IIE6 

3 

9.31186 

3 

1 

9.65736 

2 

8 65866 

2 

2 

9.31166 

2 

9.31496 

2 

3 

9.59976 

2 

9.51916 

2 

4 

9. 39616 

2 

9.49926 

2 

5 

9.52366 

2 

9.55996 

2 

6 

9.52446 

2 

9 56436 

2 

7 

9 37296 

2 

9.41216 

2 

a 

9.24336 

2 

0.27896 

2 

9 

9. 14686 

2 

9.17226 

2 

19 

9 18336 

2 

9.22546 

2 

11 

0. l'92f 

2 

9.21886 

2 

12 

9.18756 

2 

9.25326 

2 

13 

9. 1-976 

2 

9.243X 

2 

14 

9.25916 

2 

9 37326 

2 

IS 

9.1S906 

2 

9 25896 

2 

16 

9.39596 

2 

9 52996 

2 

17 

9.595-1 

2 

9.11116 

3 

18 

9. 42776 

3 

9 86566 

3 

19 

9. 19156 

4 

0 2244E 

4 

29 

9. 11S96 

J 

9.27936 

3 

21 

9 26816 

2 

9.72976 

2 

22 

9.21146 

2 

9.61416 

2 

23 

9.19216 

2 

9. 34456 

2 

24 

9. 19486 

2 

0 49146 

2 

25 

9. 76246 

1 

9.33476 

2 

26 

9.33986 

1 

9 . 45776 

2 

27 

9.55946 

1 

9.33376 

2 

28 

9 . 42656 

1 

9 . 39256 

2 

( 

U C3P1K l 

£>♦« r » 

9 

9.31186 

3 

9. 18486 

3 

1 

9.65866 

2 

9.11876 

3 

2 

9.31496 

2 

9.45l«£ 

2 

3 

9.S191E 

2 

9. 44946 

2 

4 

9 43926 

2 

9 «72lC 

2 

5 

9.55996 

2 

9.51886 

2 

6 

9 564 36 

2 

9.52236 

2 

7 

9.41216 

2 

9.41666 

2 

a 

9.27896 

2 

9.28516 

2 

9 

9.17226 

2 

9.21296 

2 

19 

9.22546 

2 

9.21956 

2 

II 

9.21886 

2 

9 22996 

2 

12 

9.25326 

2 

9.24216 

2 

13 

9.24336 

2 

9.27956 

2 

14 

9.37826 

2 

9.31446 

2 

15 

9.25996 

2 

9. 35586 

2 

1C 
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PREFACE 


This is the User Manual for the software developed for 
the evaluation of the Wiener Spectrum and Modulation Transfer 
Function of the Autoradiographic Enhancement Process. 

This report is Volume 3 of the Final Report prepared by 
ESPEE, INC. under Contract No. NAS8- 33405 for the Space 
Sciences Laboratory of the George C. Marshall Space Flight 
Center. The NASA COR for this contract is Dr. C. A. Lundquist. 
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1. INTRODUCTION 


This manual is a description of the use of the Autoradiogr- 
aphic Data Analysis software package with the Data General 
NOVA 1200 computer system. 


In the following description, the symbol ) means a carriage 
return on the tektronix. If at any time a typing mistake is 
made, the character typed may be deleted by typing the RUBOUT 
key. Any execution of a program or command may be nullified, 
by pressing the CNTRL and A keys simultaneously. 

In case of problems, call the person in charge of the 
computer system. 


2. NOVA POWER-UP PROCEDURE 


Before the Autoradiographic Data Analysis software package 
may be used, the NOVA system must be powered up and the operating 
system placed into the computer. The procedure for doing this 
is as follows: 

1. lurn power on to the NOVA computer using the front panel 

key. 

2. Turn power on to the die* controller (mounted on the DATA 
GENERAL rack ) using the front panel key. 

3. Make sure the ^ack labelled ESPEE is mounted in the disk 
drive. 

If it is not, the pack labelled ESPEE may be mounted as 
follows : 

(a) . Remove the pack inside the drive by following 

these steps: 

(i). Raise the lid of the disk drive. 

(ii). Set the clear plastic cover onto the 
disk pack. 

(iii). Turn the handle counterclockwise until 
a pop is heard. 

(iv). Raise the disk pack out and place it on 
the bottom cover. 

(v) . Lock the bottom cover. 

(b) . Replace the pack marked ESPEE inside the drive 

by following these steps: 

(i). Unlock and remove the bottom cover from 
the disk pack. 

(ii). Place the pack in the drive. 

(iii). Turn the handle clockwise until tight. 

(iv). Remove the clear plastic cover. 

(v) . Close the lid of the disk drive. 

4. Turn power on to the disk drive by pressing the on button 
on the drive. Wait for the green light on the disk drive 
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to come on before proceeding. 

5. Turn power on to the Tektronix graphics display. 

After the screen brightens, press the CLEAR button. Flip 
the LINE/LOCAL switch to LOCAL. Press the ESC key. Press 
the ; key. Flip the LINE/LOCAL switch back to LINE. This 
sets the character size for the terminal to the smallest. 

For medium and large character sizes respectively, press 
9 and 8 instead of ; in the above procedure. 

6. Turn power on to the Hard Copy Unit. The unit requires 
about 10 minutes to warm up. 

7. Set the switches on the NOVA computer panel to octal 00033. 
(This corresponds to switches 11,12,14 and 15 up; all other 
down). Raise the RESET switch momentarily. Next, raise 

the switch 0 on the computer panel up. Now ra: 1 ' momentarily 
the PROGRAM LOAD switch on the computer panel. 

This series of operations loads the operating system into 
the computer memory. The word FILENAME? will appear on the 
Tektronix screen. 

8. Type the word SYS32K ^ . When requested by the computer, 
enter the date using the format MM/DD/YY ^,and the time 
using the format HH:MM:SS ). Computer types R. 

This completes the steps required to get the computer 
powered-up. The next series of steps are required to prepare 
for using the Autoradiographic Data Analysis pachage. 
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3. TAPE MOUNTING PROCEDURE 


The source of input data for the Autoradiographic Data 
Analysis package is the Microdensitometer output tape. The 
procedure for mounting the tape is as follows; 

1. Turn power on to the magnetic tape recorder. 

2. Pull out the lever marked :OVERRIDE 

3. Mount the tape on the available slot and turn the handle 
clockwise until tight. 

4. Thread the tape as shown in the diagram on the tape recorder 
unit, and a few feet of tape on the empty reel. 

5. Press the LOAD/UNLOAD button on the recorder. 

6. Press the FORWARD button on the recorder. 

7. When the tape stops winding, press th? REMOTE button on 
the recorder. 

The tape recorder is now ready for operation. Type 
INIT MT0 ) on the Tektronix terminal. The system is now ready 
for program execution. 
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4. PROGRAM EXECUTION 


The program can be executed with four options: (1) Without 

Filters, (2) With Low Frequency Filter Only, (.3) With High 
Frequency Filter Only, and (4) With Both Filters. For program 
execution the following steps are followed: 

1. Without Filters 

1. Type DIR ESPEE ) . Computer types R. 

2. Type SP1 ) . 

The crmputer will now ask for the user-specified inputs 
tc tne program, one at a time. These inputs are: 

FILE NUMBER, NF- 
DIFFUSE DENSITY FACTOR ,DDF- 
MAXIMUM LAG FACTOR, M- 
IIEIGHT OF SCANNING SLIT.XL- 
WIDTH OF SCANNING SLIT.XA- 

3. After typing each value, type ^ . When all values are 
specified, the computer will halt after program execution 
is complete. 

2 . With Low Frequency Filters Only 

1. Type DIR ESPEE ). Computer types R. 

2 . Type SPL ) . 

The computer will now ask for user-specified inputs 
to the program, one at a time. These inputs are: 

FILE NUMBER, NF» 

DIFFUSE DENSITY FACTOR, DDF* 

MAXIMUM LAC FACTOR, M- 
HEIGHT OF SCANNING SLIT.XL- 
WIDTH OF SCANNING SLIT.XA* 

LOW FREQUENCY FILTER PARAMETER , P= 

3. After typing each value, type J . When all values are 
specified, the computer will halt after program execution 
is complete. 
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With Hl^n Frequency Filters Only 

1. Type DIR ESPEE ) . Computer types R. 

2. Type SPH ) . 

The computer will now ask for user -specif ied inputs 
to the program, one at a time. These inputs are: 

FILE NUMBER, NF" 

DIFFUSE DENSITY FACTOR, DDF- 
MAXIMUM LAG FACTOR, M- 
HEIGHT OF SCANNING SLIT.XL- 
WIDTH OF SCANNING SLIT.XA- 
HIGH FREQUENCY FILTER PARAMETER , Q- 

3. After typing each value, type ). When all values are 
specified, the computer will halt after program execution 
is complete. 


4. 


With Both Filters 

1. Type DIR ESPEE ) . Jomputer types R. 

2 . Type SP2 ) . 

The computer will now ask for user-specified inputs 
to the program, one at a time. These inputs are: 

FILE NUMBER, NF" 

DIFFUSE DENSITY FACTOR, DDF" 

MAXIMUM LAG FACTOR, M* 

HEIGHT OF SCANNING SLIT.XL- 
WIDTH OF SCANNING SLIT.XA" 

LOW FREQUENCY FILTER PARAMETER, P- 
HIGH FREQUENCY FILTER PARAMETER, Q« 

3. After typing each value, type ) . When all values are 
specified, the computer will halt after program execution 
is complete. 
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5. OBTAINING A HARD COPY 


The results of the program execution may be recorded by 
using the Hard Copy Unit. A copy of the Tektronix display 
may be obtained by pressing the AUTO PRINT/COPY switch on the 
Tektronix keyboard to COPY. 

In case the results take more than one page, automatic 
copying of the results may be accomplished as follows: 

1. Set the AUTO PRINT/ COPY switch to AUTO PRINT. 

2. Set the MARGIN CONTROL switch on the Tektronix terminal 
to 2. 
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6. SYSTEM SHUT DOWN PROCEDURE 

The procedure for shutting down the system is as follows: 

1 . Type RELEASE MT6 ) . 

2. Type RELEASE D?6 ) . This will cause the computer to halt. 

3. Turn power off to the hard copy unit. 

4. Turn power off to the Tektronix terminal. 

5. Turn power off co the disk drive by pushing the ON button. 

6. Press the LOAD/UNLOAD button on the tape recorder. Manually 

rewind the tape and remove it from the tape recorder. 

7. (Caution: Before executing this step, make sure that the 
disk has stopped spinning.) Turn power off to the disk 
controller. 

8. Turn power off to the NOVA computer. 

The shut down procedure is now complete. 


[ 
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7. PROGRAM LIMITATIONS 


The following are some of the limitations of the program. 

In actual practice, these limitations are not likely to be 

violated. 

1. The maximum no. of data points on a record on the microde- 
nsitometer tape is 1000. 

2. The maximum no. of records per file on the microdensitometer 
tape is 50. 

3. The maximum value of the Maximum Lag Factor, M is 100. 

4. For comparison with the microdensitometer printout, at 
most 300 data points will be printed from the first data 
record. 
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0.58831 

2 

0.633IE 

2 

1.78 

1.76 

1.81 

1 .85 

1.84 

1.85 

1.83 

1.81 

1.80 

1.80 

7 

0.3970C 

2 

0. 43d8£ 

2 

1.77 

1.78 

1.79 

1.77 

1.81 

1.79 

1.76 

1.79 

1.77 

1.80 

8 

0.2616E 

2 

0.29931 

2 

1.84 

1 .85 

1.87 

1 .84 

1.75 

1.79 

1.82 

1.82 

1.78 

1.76 

9 

0.1585E 

2 

0.1373E 

2 

1.78 

1.78 

1.80 

1 .80 

1.80 

1.78 

1.79 

1.81 

1.79 

1.88 

10 

0 . 1834E 

2 

0.225SE 

2 

1.77 

1 . 75 

1.79 

1.81 

1.77 

1 .77 

1.78 

1.74 

1.74 

1.77 

11 

0.1460E 

2 

0.1877E 

2 

1.79 










12 

0. 1650E 

2 

0.2228E 

2 
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0.1496E 

2 

0.2132E 

2 


9.1604E 

0.10341 
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17 

0.7324E 

1 

0 . 1 366E 
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18 

0.8074E 

1 

0. I634E 

2 

41 

2 


19 

0. 362 1 E 

1 

O.S0O2C 

1 

42 

2 


20 

0.323SE 

1 

0.7857E 

1 

.3 

2 


21 

0.2194E 

1 

0.59O0E 

1 

44 

34 


22 

0.2383E 

1 

0.8646E 

1 

45 

206 


23 

0.17661 

1 

0.5960E 

1 

46 

349 


24 

0.1882C 

1 

0.7 206E 

1 

47 

206 


2S 

0.1269E 

1 

0.S572E 

1 

48 

41 


26 

0.181 SE 

1 

0 . 92 3 IE 

1 

49 

1 


27 

0.U60E 

1 

0.69171 

1 

RAU MEAN • 

0. 14S498E 

4 

28 

0. 1087E 

1 

8.7709E 

1 

MEAN OF NORMALIZED DATA- 

0 . 1 27000E 1 

K 

GKC3PIK 

• 

CHiKI 

MEAN SQUARE 

VALUE • 0 . 

•687189E -3 

0 

0.3I33E 

3 

0.19161 

3 


r»ni«\iajAi<urg<unjfviajftjf\jnjnjajf\jnjai 


iS53l8*SSSSlSS=5SS*5)«;*?5a5?a?SS 

•«l/»YlAU)l/tVI*t(\IAI(M(UM 


< r- r* h- • A U) <i) r-r» r- 


AiniiUfUf\ininiiunii\ji\i(\iivj(UMrg<\i«u 



MninTtfi^hvot^^funtt/KAhoiQif 

— — — — — — — — ^Aj(\jfO<UnjOJ«VAJf\J 


III-ll 


STOP 


% 


\ 


SPL 

FILE NUNBER.NF-S 
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HEIGHT OF SCANNING SLIT.XL-IO00 
WIDTH OF SCANNING SLIT.XA-25 
LOU FREQUENCY FILTEP PaRaNETER.P- 30 
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NUMBER Of FILES SUPPED- 4 
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NUMBER OF SAMPLES- 281 
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V COORD. NATE (MICRONS)- 4 

DELTA X iHI (MICRONS*- 18 

DATA IS BEING FORMATTED 
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1.68 

1.69 

1.80 

1.83 

1.82 

1.81 

1.76 

1.79 

1.85 

1.82 

1.84 

1.87 

1.87 

1.81 

1.77 

1.79 

1.83 

1.82 

1.82 

1.84 

1.86 

1.82 

1.81 

1.84 

1.82 

1.84 

1.81 

1.86 

1.89 

1.89 

1 .85 

1.82 

1.80 

1.79 

1.82 
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1.81 
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1.79 

1.82 

1.81 

1.84 
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1.88 
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1.79 

1.79 
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1.89 
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1.81 
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1.84 
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1.79 
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1.81 
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1.79 

1.79 
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1.81 

1.79 

1.73 

1 .80 

1.82 

1 .85 
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1.78 

1.76 

1.74 

1 . 75 
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1.79 

1.78 

1.76 

1.77 

1.79 

1.82 

1.78 

1.78 

1.78 

1.78 

1.88 
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1.80 
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1.77 
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1.82 

1 .86 
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1.81 

1.78 

1.82 
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1.81 

1.80 

1.81 
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1.78 

1.79 
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1.76 
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1.82 

1.82 
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1.81 
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1.81 

1 .82 
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1.83 
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1.77 

1.81 

1.77 
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1.78 

1.79 
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1.79 

1.78 

1.77 
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1.78 

1.76 
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1.77 

1.78 

1.79 

1.77 
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1.79 
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1 .87 
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NORM. R-MAT 




0 

0.4714E 

-3 

0.1000E 

1 




1 

0.2918E 

-3 

0.6189E 

• 




2 

0.97306 


0.2064E 

0 




3 

0.2736E 

-4 

0.5803E 

-f 




4 

-0.1S73E 

-4 

-0.3336E 

• 1 




5 

-0. 4244E 

• 4 

-0.9001E 

-f 




6 

-0.4077E 

-4 

-0.8648E 

-1 




7 

-0.2055E 

-4 

-0 .43586 

- J 




8 

0.4914E 

-5 

0.1042E 

- J 




9 

0 . 224 IE 

-4 

0.47S4E 

-( 




10 

0.9430E 

-5 

0.2000E 

-J 




11 

-0.1482E 

-4 

-0.3144C 





12 

-0.1328E 

-4 

-0.2816E 


1.78 

1 .82 

1.79 

13 

-0.I152E 

-4 

-0.24436 

-1 

1.80 

1 .79 

1.81 

14 

0.2003E 

-/ 

-0.4249E 

- 1 

1.81 

1.80 

1.85 

15 

-0.1S43C 

-4 

-0.3273E 

-1 

1.78 

1.81 

1.85 

16 

-0.4060C 

-S 

-0 8612E 

-2 

1.79 

1.80 

1.81 

17 

-0.2292E 

-6 

-0.4862E 

-3 

1 83 

1.85 

1.82 

18 

0.8791E 

-5 

-0.1780E 

-| 

1.81 

1 .82 

1.85 

19 

-0.20726 

-4 

-0. 4336E 

-1 

1.30 

1 .81 

1.86 

20 

-0.2S93E 

-4 

-0.5500E 


1.85 

1.83 

1.84 

21 

-0.24S6E 

-4 

-0.5209E 

-1 

1.79 

1.78 

1.81 

22 

-3.2420E 

-4 

-0.51 326 

-I 

< .80 

1.79 

1.80 

23 

-0.3434E 

-4 

-0.728SE 


1 .77 

1.79 

1 .78 

24 

-0.4920E 

-4 

-0.10446 

0 

1.74 

1.76 

1.77 

2S 

-0.5663E 

-4 

•0.12016 

0 

1.77 

1.77 

1.80 

26 

-0.7060E 

-4 

-0.1 4986 

0 

1.78 

1.84 

1.83 

27 

-O.7054E 

-4 

-0.14966 

0 

1.77 

1 .80 

1.77 

28 

0.4194E 

-4 

-0.8897E 

-1 

1.82 

1 .80 

1.81 
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GKC** 1 " « > 


1.81 

1.82 

1.80 

0 

0.6046E 

! 

«• 60 4 fct 

1 

1.82 

1.86 

1.86 

t 

0.6598E 

a 

2 .-SUE 

2 

1.86 

1 .88 

1.88 

2 

0.3026E 

a 

« ;iS»f 

2 

1.85 

1 .83 

1.79 

3 

0.S156E 

a 

o.:...~ 

2 

1.80 

1.81 

1.80 

4 

0.39S5C 

a 

0.4086E 

2 

1.77 

1.77 

1.80 

S 

0.45196 

a 

0.4754E 

2 

i.at 

1 .80 

1.80 

6 

0. 4379E 

a 

0.4712E 

2 

1.79 

1.77 

1.80 

7 

0.3648E 

a 

A.4033E 

2 

1.82 

1 .78 

1.76 

8 

0.2371E 

a 

O.2703E 

2 

1.81 

1 .79 

1 .80 

9 

0.1355C 

a 

0.16OOE 

2 

1.74 

1.74 

1.77 

10 

0.17476 

a 

9.2148E 

2 




11 

0. 1095E 

a 

0.1407E 

2 




12 

0. 1540E 

a 

0.2030E 

2 




13 

0.1479E 

a 

0.21086 

2 




14 

0. 12ISE 

a 

0. 1?37E 

2 




IS 

0.6I69E 

i 

. 995 1 E 

1 


TOTAL NO. OF 

CLASSES- 100 

16 

0.8443E 

1 

0.I460C 

2 

CLA .S NO. 

NO. OF DATA POINTS 

17 

0.6477E 

1 

0. 1238E 

2 

«! 

2 

18 

0.8799E 

1 

0.I781E 

2 

42 

2 

19 

0.34O4E 

1 

0.7S2IE 

1 

43 

2 

20 

9.4602E 

1 

0.1II8E 

2 

44 

34 

21 

0. 1429E 

1 

0.3843E 

1 

45 

206 

22 

0.3939E 

1 

0.U8IE 

2 

46 

349 

23 

0.14O9E 

1 

0.475SE 

1 

47 

206 

24 

0.2528E 

1 

0.9683E 

1 

48 

41 

£S 

0.7372E 

0 

0. 32366 

1 

49 

1 

26 

0.1925E 

1 

0.9791E 

1 

RAU MEAN • 

0.14S498E 4 

27 

0.620)6 

0 

0.370OC 

1 

MEAN OF NORMALIZED DATA* 0.127000E 1 

28 

0.1759E 

1 

0.1247E 

2 
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i 

bO 

1 
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1 

89 

1 

83 

1 

.82 


.81 

1 

.78 

1 
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1 

.79 
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i 

?9 

I 
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1 
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t 

84 

1 
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i 
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1 

.89 

1 
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1 
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1 
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1 

82 

1 
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i 

84 

1 
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1 

.89 

1 
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• 
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1 

81 

1 
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t 

82 

1 
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1 

.81 

1 
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1 

.81 

1 
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t 
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89 
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.89 
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i 
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1 

82 

t 

89 

1 

81 

t 

.82 

i 

82 

t 
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1 
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1 
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i 

82 

1 
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1 
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1 
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1 

84 

l 

.89 

1 

81 

1 
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1 
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1 
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t 

84 

1 

79 

t 

79 

i 

.82 

t 

89 

1 

81 

1 
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I 

.85 

1 

89 
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t 
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1 
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1 

85 

1 
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t 
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1 

78 

i 

78 

1 

.77 

1 

.77 

t 

.89 


>0 

• 

83 

1 

83 

t 

89 

1 

89 

1 

81 

i 

.77 

1 

.78 

1 

84 

1 

.83 


s; 

• 

82 

1 

85 

t 

83 

1 

84 

1 

81 

i 

78 

1 

.77 

1 

89 

1 

.77 



t 

87 

I 

83 

1 

81 

1 

89 

1 

81 

l 

82 

1 

82 

1 

86 

1 

.81 


5 ! 

i 

75 

1 

73 

t 

78 

1 

79 

I 

75 

l 

76 

1 

31 

1 

82 

1 

80 


*# 

i 

82 

1 

82 

l 

84 


83 

I 

89 

i 

81 

1 

82 

1 

86 

1 

86 


* * 


85 

1 

84 

1 

84 

1 

88 

t 

88 

i 

84 

1 

86 

1 

88 

1 

88 


1 » 

• 

8 1 

] 

82 

1 

84 

1 

84 

1 

83 

l 

83 

1 

85 

1 

83 

1 

.79 


• » 


77 

t 

81 

1 

.77 

1 

84 

1 

.88 

l 

.81 

1 

89 

1 

81 

1 

89 




79 

1 

83 

1 

79 

I 

73 

1 

.77 

t 

76 

1 

.77 

1 

.77 

1 

89 


rs 


76 

I 

81 

I 

8S 

1 

84 

I 

.85 

i 

83 

1 

.31 

1 

80 

1 

39 


• .• 


78 

1 

79 

1 

.77 

1 

81 

t 

.79 

t 
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41 

2 


2 

4) 

a 

44 

34 

4S 

296 

4ft 

349 

4 t 

296 

44 

41 

44 

1 
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DATA 
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ROOT MEAN SQUARE UALUE • 9.22943SE 

AUTOCORRELATION IS BEING COMPUTED 

-1 

R 

R-MAT 

NORM. R-MAT 

• 

9.486SC -3 

9.1909C 


1 

9.4I73E -3 

9.SS78E 


2 

9.39S9E -3 

9.6269E 


3 

9.2919E -3 

9.41 SRC 


4 

9.1415E -3 

9.2999E 


S 

9.11S7E -3 

9.2377C 


6 

9.I138E -3 

9.2349C 


7 

9.1395E -3 

9.2683C 


8 

9.1524C -3 

9.3132E 


9 

9.1668E -3 

9.3429E 


19 

9.16S7E -3 

9.3497C 


11 

9.1SS8E -3 

9.32C3E 


12 

9.14241 -3 

9.2927C 


13 

9.1336E -3 

9.274SC 


14 

9.12951 -3 

9.2661E 


IS 

9.1329E -3 

9.2733E 


16 

9.I376C -3 

9.2829C 


17 

9.1399E -3 

9.28761 


18 

9.1387E -3 

9 . 285 IE 


19 

9. 138 IE -3 

9 . 28391 


29 

9.13S7E -3 

9.2789E 


21 

8.1393E -3 

9.2677E 


22 

0.1224E -3 

9.25161 


23 

B.U38E -3 

9.2349E 


24 

9.1941E -3 

9.21391 


2S 

9.9SS6E -4 

9.1964E 


26 

0.8979E -4 

9.1846E 


27 

0.89S3E -4 

9.18491 


28 

9.9297E -4 

9. 191 IE 



THE DENOMINATOR OF sine function ielou I.H-M ENCOUNTERED in INF 
THE DENOMINATOR OF SINE FUNCTION BELOW t.OC-96 ENCOUNTERED IN INF 
I HE DENOMINATOR Of SINE FUNCTION BELOW 1.9C-96 ENCOUNTERED IN INF 
THE DENOMINATOR OF SINE FUNCTION BELOW I.OE-96 ENCOUNTERED IN INF 
THE DENOMINATOR Of SINE FUNCTION BELOW I.OE-96 ENCOUNTERED IN INF 
Ttf DENOMINATOR OF SINE FUNCTION BELOW I.H-W ENCOUNTERED IN INF 
THE DENOMINATOR OF SINE FUNCTION BELOW I.OE-06 ENCOUNTtRED IN INF 
TIC DENOMINATOR OF SINE FUNCTION BELOW l.OE-96 ENCOUNTERED IN INF 
THE DENOMINATOR OF SINE FUNCTION BELOW 1.9C-96 ENCOUNTERED IN I»C 
THE DENOMINATOR OF sine FUNCTION BELOW I.OE-96 ENCOUNTERED IN INF 
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GKC Pi K 1 


9 

9.3113E 

3 

9. 31 18E 

3 

1 

0.6S18E 

a 

9.6S73C 

a 

2 

9.39121 

2 

9.311SE 

2 

3 

9.4723E 

a 

9.S99~E 

2 

4 

9.345SE 

2 

9. 396 1C 

2 

5 

9.4229E 

o 

9.S236E 

2 

6 

9.3829E 

1 

9.S244E 

2 

7 

9.241SC 

2 

0.3729E 

2 

8 

9. I368C 

a 

9.2438E 

a 

• 

9 . 69921 

t 

9.14581 

2 

1* 

9.7I97C 

i 

9.18331 

2 

11 

9.S214E 

t 

9.1792E 

2 

12 

9.4JS9E 

i 

0.187SE 

a 

13 

9. 23411 

l 
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PREFACE 


This is the Description and Program Listings of the 
modification of Data Analysis Software developed for the 
evaluation of the Wiener Spectrum and Modulation Transfer 
Function of the Autoradiographic Enchancement Process. 

This report is an addendum to th« three-volume report 
No. 79-01 of Contract No. NAS8-33^05 and is the Final Report 
prepared by ESPEE, INC. under Contract No. P.0. H-30573B for 
the Space Sciences Laboratory of the George C. Marshall 
Space Flight Center. 

The NASA COR for this Contract is Mr. Ray Hembree. 
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1. MODIFICATION OF DATA ANALYSIS SOFTWARE 


The following modifications have been done 
to data analysis software developed under an earlier contract 
No. NAS8- 33405 : 

I. 1 Individual Record Processing 

The data from each record in a file on micordensitometer 
tape is individually processed. All processes, namely compute 
histogram, compute average, normalize data .calculate zero-mean 
data, filter data for low frequency, filter data for high frequence , 
compute standard deviation compute autocorrelation, compute 
power spectral density function, apply inverse filter transform, 
apply aperture compensation and smoothe power spectral densities 
are applied co each record separetely. These processes were 
formerly applied to all records in a file combine-* into one 
record. 


1.2. Averaging Results of all Records 

After the autocorrelation and smoothed power spectral 
density function arrays are computed for each record, an 
average of these arrav.3 is computed. 


1.3. Computing Standard Deviation and Error 


The standard deviation and error are computed for the 
averaged autocorrelation and smoothed power spectral densitv 
arrays. The standard deviation is computed as follows 


Standard Deviation 


r ^ 1 d 2 )/M] ;i 
L nh n 


where M is the max lag factor and is the 

The error is computed 
Error ■ Standard Deviation ' (N)" 5 
where N is the Mo. of records. 


M-th data point 
as follows: 
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1.4. 


Four Cases 


The above modifications have been applied to all four 
cases, namely without filters, with low frequency filter 
only, with high frequency filter only and with both filters. 
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2. SOFTWARE MODULES AFFECTED 


The following software modules were changed: 

SP1 

SP2 

SPL 

SPH 

A new module named LAST has been added to the software. 
This nodule imputes r he standard deviation and error for the 
averaged autocorrelation and power spectral density functions. 

This module is the last .iodule called by the main nodules SPL, 

SP2, SPL and SPH. 

The modifications have been made in a manner such 
that the use of the program is unchanged from the way described 
in the users manual. 


3. PROGRAM LISTINGS 


The program listings for the modified modules 
and the new module are given in the following pages. Also 
included are the updated list of modules and code for loading 

compiled modules. 


TVPE SPLIST 

FILE SPLIST 

SP1A-RLDR/H FOR SP1 

SP2A-RLDR/H FOR SP2 

SPLA-RLDR/H FOR SPL 

SPHA-RLDR/fl FOR SPH 

SP23A-RLDR/H FOR SP23 

SP24A-RLDR/H FOR SP24 

SPl-flAlH UITHOUT FI17ERS 

SP2-HAIN WITH FILTERS 

SPL-HAIN WITH LOU FREQ. FILTER 

SPH-P1AIN UITH HIGH FREQ. FILTER 

SP3-PTDT 

SP4-SICPFL 

SP5-RDFL 

SP6-HDR 

SP7-C0RD 

SP8-RFHT 

SP9-GBVTE 

SPie-HEANS 

SP11-SIGMA 

SP12-AUT0C0R 

SP13* ( DELETED ) 

SP14-HIST 

SP15-DCLS 

SP16-AUG 

SP17-NRMLZ 

SP18-ZHEAN 

SP19-PSDF 

SP20-SPDF 

SP21-INF 

SP22-AC 

SP23-LFF 

SP24-HFF 

SP2S-AC1 

SP26-LAST 


TYPE SPl A 

RLDR/P1 SP1 SP3 SP4 SP5 SPS SP7 SP8 *P9 SP10 SPU SP12 SP14 SP15 ~ 
SP16 SP17 SP18 SP19 SP25 SP20 SP26 FORT. LB 
R 

TYPE SP2A 

RLDR/fl SP2 SP3 SPA SP5 SPG SP7 SP8 SP9 SP10 SP11 SP12 SP14 SP1S ~ 
SPl 6 SP17 SP18 SP1 9 SP21 SP22 SP20 SP26 FORT. LB 
R 

TYPE SPLA 

RLDR/fi SPL SP3 SPA SPS SPG SP7 SPS SP9 SP10 SP11 SP12 SP14 SP15 ~ 
SPiG SP17 SP18 SP19 SP21 SP25 SP20 SP26 FORT. LB 
R 

TYP" SPHA 

RLDR/H SPH SP3 SP4 SP5 SPG SP7 SP8 SP9 SP10 SPli SP12 SP14 SP15 * 
SPl 6 SP17 SPl 8 SP19 SP21 SP22 SP20 SP26 FORT. LB 
R 

TYPE SP23A 

RLDR/H SP23 FORT. LB 

R 

> TYPE SP24A 

i RLDR/H SP24 FORT. LB 

R 


o o ooooooooooooooo 


TYPE SPl 

FILE SP1-MAIN WITHOUT FILTERS 
PROGRAMMER NAME* S.P. SINGH 
REUISIONs 3, MULTIPLE RECORDS/FILE 
DATE 17/31/ 79 

PURPOSE * THIS MODULE ACCEPTS USER INPUTS AND CALLS ALL 
SUBSEQUENT MODULES 

RRH IS AUTOCORRELATION UALUES OF DATA 
GKC IS POWER SPECTRAL DENSITY FUNCTION 
GKH IS SMOOTHED POWER SPECTRAL DENSITY 
GKCP IS THE INUERSE FILTER FUNCTION 
GKC3P IS THE APERTURE COMPENSATION FUNCTION 
ITEMP.XTEMP ARE BUFFERS 
NR IS THE NUMBER OF RECORDS IN THE FILE 
NW IS THE NUMBER OF WORDS PER RECORD 
NP IS THE NUMBER OF POINTS PER RECORD 
INTEGER P.Q 

DIMENSION RRH< 101 >,ITEMP< 1000), 

1 GKC( 101 ),GKH( 101 ) 

DIMENSION GICCP ( 101 ),GKC3P( 101 ) 

DIMENSION NU(50 ),NP(50 ) 

DIMENSION XTEMP( 1000 ) 

DIMENSION RR( 101 ),GH( 101 ) 

ACCEPT 'FILE NUMBER. NF- ■ ,NF 
ACCEPT 'DIFFUSE DENSITY FACTOR. DDF-', DDF 
ACCEPT 'MAXIMUM LAG FACTOR, M-'.M 
ACCEPT 'HEIGHT OF SCANNING SLIT, XL-', XL 
ACCEPT 'WIDTH OF SCANNING SLIT.XA-'.XA 
ACCEPT 'LOW FREQUENCY FILTER PARAMETER, P-' ,P 
ACCEPT 'HIGH FREQUENCY FILTER PARAMETER, Q-',Q 
P-0 
Q-0 

CALL PTDT (ITEMP.NF.H.NR.NU.NP) 

C TRANSFER RAW DATA RECORDS INTO A FILE SRD 

CALL FOPEN( 3, 'RD' ,2000 ) 

CALL FOPEN(2,'SRD',2000> 

DO 10 I - 1 , NR 

CALL READR(3, I , ITEMP, 1 , IER ) 


CALL URITR<2,I,ITEHP,1,IER) 

10 CONTINUE 

CALL FCLOS ( 2 ) 

CALL FCLOSO) 

NRK-i 

C CALL FOPENM, 'PAR* ) 

C WRITE EINARY (4) NP,P # G»NRIC 

C CALL FCI.OS( 4 ) 

A-0 

D-320C 

NC-lvO 

DO 15 J-i,101 
RR( J )■©. 

GH( J)-0. 

15 CONTINUE 

DO 20 I ■ 1 # NR 
TYPE * * 

TYPE 'RECORD NO.- # I 
TVPE * * 

CALL HIST ( ITEMP, A, B,NC,NRK,NP ) 

CALL MEANS ( ITENP,XTEI1P,NP,r*PK # DDF ) 

C CALL SUAP< -SP23.su* , IER ); CALL LFF(XTEHP # TEMP2,TEHP3,NP # P,NRIC ) 

C CALL SUAP( -SP24.SU- , IER )j CALL HFF<XTEMP,TEf1P2,TEnP3,NP,Q,NRIC > 

CALL SIGNA(XTENP,NRIC,NP) 

CALL AUTOCOR ( XTEMP , NP, NRK , 11 , RRH, P , Q ) 

HP1-M+1 

DO 30 J-1,MP1 

RR ( J ) "RRH ( J )+RR ( J ) 

30 CONTINUE 

CALL P5DF(RRH,H,H,XL # GICC ) 

C CALL INF(GKC,GICCP,0,n) 

CALL AC1 (GKC # GKC3P,H,H,XA ) 

CALL SPDF ( f1,GKC3P, GKH ) 

DO 40 J-1.NP1 
GH( J )-GH( J)+GKH( J ) 

40 CONTINUE 

CALL DELETECRD- ) 

CALL DELETE ( -NRb* ) 



> 

I 


20 

50 

100 

60 

70 

R 


CALL DELETE CZRD* > 

CALL DELETE < *ZP* > 

CALL DELETE < *TD* ) 

CALL FOPEH(3.*RD',2000) 

CALL F0PEH(2, *SRD #2000 > 

IF (I.GE.NR)GO TO 20 

CALL + READR(2,Il,ITEhP.l.IER> 

CALL URITR(3.1.ITEHP,1,IER> 

CALL FCL0SC2 ) 

CALL FCL0S(3) 

CONTINUE 

CALL DELETE C 'PAR’ ) 

DO 50 J-1,NP1 
RR<J>-RR<J>'M* 

GHCJ)-GH(J)/NR 

TYPE' N R E AVERAGE R-HAT" 

DO 60 J-1,HP1 

IR-J-1 

WRITE (10,100)IR#RR ( J J 
FORMAT ( 1X#I3#8X,E11.4> 

t?pe* NUE average GH(K)* 

DO 70 J-l,f1Pl 

IR-J-1 

URITE( 10# 100 > IR,GH(J) 

CONTINUE 

CALL LAST <RR,GH,HP1#NR ) 

STOP 

END 



noooooooooooooo 




TYPE SPL 

FILE SPL-MAIN UITH LOU FREQ. FILTER 
PROGRAMMER NAME * S.P. SINGH 
REVISION t 3 , MULTIPLE RECORDS/FILE 
DATE 1 8/1-4/79 

PURPOSE t THIS MODULE ACCEPTS USER INPUTS AND CALLS ALL 
SUBSEQUENT MODULES 

RRH IS AUTOCORRELATION VALUES OF DATA 
GKC IS POWER SPECTRAL DENSITY FUNCTION 
GKH IS SMOOTHED POUER SPECTRAL DENSITY 
GKCP IS THE INVERSE FILTER FUNCTION 
GKC3P IS THE APERTURE COMPENSATION FUNCTION 
ITEMP,XTEHP ARE BUFFERS 
NR IS THE NUMBER OF RECORDS IN THE FILE 
NU IS THE NUMBER OF UORDS PER RECORD 
NP IS THE NUMBER OF POINTS PER RECORD 
INTEGER P,Q 

DIMENSION RRH(101 >, ITEMP( 1000 >, 

1 GKC( 101 ),GKH( 101 ) 

DIMENSION GKCP( 101 ),GKC3P(101 ) 

DIMENSION NU(50 ),NP(50 ) 

DIMENSION XTEMP( 1000 ) 

DIMENSION RR( 101 ),GH( 101 ) 

ACCEPT "FILE NUMBER, NF- • ,NF 
ACCEPT 'DIFFUSE DENSITY FACTOR, DDF- ', DDF 
ACCEPT 'MAXIMUM LAG FACTOR, M-'.H 
ACCEPT 'HEIGHT OF SCANNING SLIT, XL-', XL 
ACCEPT 'UIDTH OF SCANNING SLIT,XA-',XA 
ACCEPT 'LOU FREQUENCY FILTER PARAMETER, P- ' ,P 
C ACCEPT 'HIGH FREQUENCY FILTER PARAMETER, Q-',Q 

Q-0 

CALL PTDT ( ITEMP,NF,H,NR,NU,NP ) 

C TRANSFER RAU DATA RECORDS INTO A FILE SRD 

CALL FOPENO, 'RD',2000) 

IS-1 

CALL FOPEN< IS, 'SRD ',2000 ) 

DO 10 I-1,NR 

CALL READR( 3, I , ITEMP, 1 , IER ) 


10 


15 


> 

t 


C 


C 


30 

C 


CALL URITR(IS,I, ITENP, 1 , IER ) 

CONTINUE 
CALL FCLOS(IS) 

CALL FCLOSO) 

NRIC-1 

CALL F0PEN<4, 'PAR' ) 
uriTE BINARY (4) NP»P,Q,NRK 
CALL FCL0SC4) 

A-0 
B-3200 
NC* 100 

DO 15 J-1,101 
RR(J)-0. 

GH(J)-0. 

CONTINUE 
DO 20 I • 1 , NR 
TYPE ' ' 

TYPE 'RECORD NO.', I 
TYPE ' ' 

CALL HIST ( ITENP, A,B,NC,NRK,NP ) 

CALL MEANS C ITEMP,XTENP,NP,NRK, DDF ) 

CALL SUAP('SP23.SU',IER); CALL LFF C XTEHP, TENP2, TEMP3, NP. P, NRK ) 
TYPE ' IER* ' , IER 
CALL DELETE CZRD') 

CALL RENAN ( 'ZP' , 'ZRD* , IER ) 

CALL SUAP ( ' SP24 . SU' , IER ) j CALL HFF(XTENP,TENP2,TENP3,NP,Q,NRK ) 
CALL SIGNA<XTENP,NRK,NP) 

CALL AUTOCOR ( XTENP , NP , NRK , N , RRH , P , Q ) 

NP1-H+1 

DO 30 J-1,NP1 

RR( J)-RRH( J)+RR( J ) 

CONTINUE 

CALL PSDF <RRH,M,H,XL,GKC ) 

CALL INF ( GKC, GKCP , Q , H ) 

CALL AC1 ( GKC, GKC3P, N, H, XA ) 

CALL SPDF (N, GKC3P, GKH ) 

DO 40 J* 1 ,NP1 
GH(J)-GH(J)+GKH(J) 


40 CONTINUE 

CALL DELETE( 'RD' ) 

CALL DELETECNRD' ) 

CALL DELETE ( 'ZRD* ) 

CALL DELETECZP' ) 

CALL DELETE('TD') 

CALL F0PEN(3, *RD',E000) 

CALL FOPENUS, 'SRC*, 2000) 

IF (I.GE.NRJGO TO 20 
Il-I+l 

CALL RE ADR ( IS, II , I TEMP, 1 , IER ) 
CALL URITR(3,1,ITEF.P,1,IER) 

CALL FCLOS(IS) 

CALL FCLOSO) 

20 CONTINUE 

CALL DELETE ( 'PAR 1 ) 

DO 50 J-1,HP1 
RR( J)-RR(J)/NR 
GH( J)-GH(J)/NR 
50 CONTINUE 

TYPE' R AVERAGE R-HAT 

DO 60 J-1,MP1 

IR-J-1 

URITE (10/100 )IR,RR( J ) 

100 FORMAT ( IX, 13, 8X, Ell. 4) 

6* CONTINUE 

TYPE' K AVERAGE GH(K) 

DO 70 J-1,MP1 

IR-J-1 

URITE ( 10, 100 ) IR, GH( J ) 

70 CONTINUE 

CALL LAST (RR,GH,MP1 ,NR ) 

STOP 

END 

R 


ooooooooooooooo 


TYPE SPH 

FILE SPH-flAIN UITH HIGH FREG. FILTER ONLY 
PROGRAMMER NANEt S.P. SINGH 
REUISI0N*3, MULTIPLE RECORDS/F ILC 
DATE *8/20/79 

PURPOSE* THIS MODULE ACCEPTS USER INPUTS AND CALLS ALL 
SUBSEQUENT NODULES 

RRH IS AUTOCORRELATION UALUES OF DATA 
GKC IS POWER SPECTRAL DENSITY FUNCTION 
GKH IS SMOOTHED POWER SPECTRAL DENSITY 
GICCP IS THE INVERSE FILTER FUNCTION 
GKC3P IS THE APERTURE COMPENSATION FUNCTION 
ITEMP,XTEHP ARE BUFFERS 
NR IS THE NUMBER OF RECORDS IN THE FILE 
NW IS THE NUMBER OF WORDS PEP RECORD 
NP IS THE NUMBER OF POINTS PER RECORD 
INTEGER P,Q 

DIMENSION RRH( 101 ), I TEMP ( 1000 ), 

1 GKC( 101 ),GKH( 101 ) 

DIMENSION GtCCP( 101 ),GKC3P( 101 ) 

DIMENSION NU( 50 ) , NP ( 50 ) 

DIMENSION XTEMPU000) 

DIMENSION RR( i0i ),GH( 101 ) 

ACCEPT 'FILE NUMBER, NF- • ,NF 
ACCEPT 'DIFFUSE DENSITY FACTOR, DDF ', DDF 
ACCEPT 'MAXIMUM LAG FACTOR, M-',M 
ACCEPT 'HEIGHT OF SCANNING SLIT, XL-', XL 
ACCEPT 'WIDTH OF SCANNING SLIT,XA-',XA 
C ACCEPT 'LOW FREQUENCY FILTER PARAMETER, P- * ,P 

P-0 

ACCEPT 'HIGH FREQUENCY FILTER PARAMETERS- ' ,Q 
CALL PTDT (ITEMP,NF.H,NR,NW,NP) 

C TRANSFER RAU DATA RECORDS INTO A FILE SRD 

CALL F0PEN(3, 'RD',2000) 

IS-1 

CALL FOPENdS, 'SRD', 2000) 

DO 10 I • 1 , NR 

CALL READR(3,I, ITEMP.l, IER ) 


-14 


ie 


15 


> 


c 


30 


CALL UR I TR (IS, I, ITEMP, 1 , IER ) 

CONTINUE 
CALL FCLOS(iS) 

CALL FCLOSO) 

NRK-1 

CALL F0PEN(4, 'PAR* ) 

WRITE BINARY (4) NP,P,<2,NRK 
CALL FCLOS( 4 ) 

A'C 

B-3200 
NO- 100 

DO 15 J-1,101 
RR(J)-9. 

GK(J)“C. 

CONTINUE 
DO 20 1*1, NR 
TYPE ' ' 

TYPE 'RECORD NO.', I 
TYPE ' • 

CALL HIST ( ITEMP, A, B,NC , NRK, NP ) 

CALL MEANS ( I TEMP , XTEMP, NP , NRK , DDF ) 

CALL SUAPC *SP23.SU*,IER)j CALL LFF( XTEMP, TENP2,TEMP3,NP,P, NRK ) 
rfli 1 PPNAM ( •7RD* a 7P* TPP 1 

CALL SUAP(*SP24.sO*,IER)j CALL HFF( XTEMP, TEMP2,TEMP3,NP,0, NRK ) 

CALL DELETE CZRD'i 

CALL RENAM ( *ZPP* , *ZRD* , IER ) 

CALL SIGHACXTEMP,NRK,NP> 

CALL AUT0C0R( XTEMP, NP, NRK, r,,RRH,P,Q) 

MP1-M+1 

DO 30 J-1,MP1 

RR( J)-RRH(J)+RR( J) 

CONTINUE 

CALL PSDF(RRH,M,H,XL,GKC ) 

CALL INF(GKC,GKCP,Q,M) 

CALL AC( GKCP,GKC3P, M, H, XA ) 

CALL SPDF ( H,GKC3P, GKH ) 

DO 40 J-1.MP1 
GH( J)-GH(J)+GKH(J) 


-15 


40 CONTINUE 

CALL DELETE! "RD* ) 

CALL DELETE! 'NRD* > 

CALL DELETE! *ZRD' ) 

CALL DELETE ! # ZP'> 

CALL DELETE! , TD* ) 

CALL FOPEN!3,'RD',2000> 

CALL FOPEN! IS, 'SRD',2000) 

IF !I.GE.NR)GO TO 20 
11 - 1*1 

CALL RE ADR! IS, II, ITENP, 1, IER ) 
CALL URITR! 3, 1 , ITENP, 1 , IER ) 

CALL FCLOS! IS ) 

CALL FCLOS! 3) 

20 CONTINUE 

CALL DELETE! 'PAR* ) 

DO 50 J-1,NP1 
RR! J)-RR! J)/NR 
GH( J )-GH( J )/NR 
50 CONTINUE 

TYPE* R AVERAGE R-HAT* 

DO 60 J- 1 , NP1 

IR-.l-l 

URITE ! 10,100)IR,RR(J ) 

100 FORMAT! IX, 13, 8X, Ell *4) 

60 CONTINUE 

TYPE* K AVERAGE GH!*)* 

DO 70 J-1,HP1 

IR-J-1 

URITE! 10, 100 ) IR,GH! J ) 

70 CONTINUE 

CALL LAST !RR,GH,NP1,NR ) 

STOP 

END 

R 
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TYPE SP2 

FILE 5P2-MAIN WITH FILTERS 
PROGRAMMER NAME J S.P. SINGH 
REUISI0NI3, MULTIPLE RECORDS/FILE 
DATE * 8/24/79 

PURPOSE * THIS NODULE ACCEPTS USER INPUTS AND CALLS ALL 
SUBSEQUENT NODULES 

RRK IS AUTOCORRELATION VALUES OF DATA 
GKC IS POWER SPECTRAL DENSITY FUNCTION 
GKH TS SMOOTHED POUER SPECTRAL DENSITY 
GKCP IS THE INVERSE FILTER FUNCTION 
G\C3P IS THE APERTURE COMPENSATION FUNCTION 
IT£RP # XTEKP ARE BUFFERS 
NR IS THE NUMBER OF RECORDS IN THE FILE 
NU IS THE NUMBER OF UOFDS PER RECORD 
NP IS THE NUMBER OF POINTS PER RECORD 
INTEGER P,G 

DIMENSION RRH(101 ), ITENP< 1000 ), 

1 GKC( 101 ),GKH( 101 ) 

DIMENSION GKCP( 101),GKC3P(10i) 

DIMENSION NU<S0).NP(50) 

DIMENSION XTEMP( 1000 ) 

DIMENSION RR( 101), GH( 101 ) 

ACCEPT 'FILE NUMBER , NF NF 
ACCEPT * DIFFUSE DENSITY FACTOR, DDF- ', DDF 
ACCEPT 'MAXIMUM LAG FACTOR, M-',M 
ACCEPT 'HEIGHT OF SCANNING SLIT, XL-', XL 
ACCEPT 'U1DTH OF SCANNING SLIT,XA-',XA 
ACCEPT 'LOU FREQUENCY FILTER PARAMETER, P- ' ,P 
ACCEPT 'HIGH FREQUENCY FILTER PARAMETER, Q- Q 
CALL PTDT ( ITEHP, NF,H, NR, NU,NP ) 

C TRANSFER RAU DATA RECORDS INTO A FILE SRD 

CALL FOPEN(3,'RD',2000) 

IS-1 

CALL FOPEN(IS, 'SRD', 2000) 

DO 10 I-t.NR 

CALL RE ADR ( 3, 1 , 1 TEMP, 1 , IER ) 

CALL URITR( IS, I , ITEMP, 1 , IER ) 


10 CONTINUE 

CALL FCLOS(IS) 

CALL FCLOSO) 

NRK-1 

CALL F0PEN(4, 'PAR’ ) 

WRITE BINARY (4) NP,P,Q.NRK 
CALL FCLOS (4 ) 

A-0 

B-3200 

NC-iee 

DO 15 J-1,101 
RR(J>-0. 

GH(J)-0. 

15 CONTINUE 

DO 20 I-1,NR 
TYPE ' ’ 

TYPE 'RECORD NO.', I 
TYPE * ' 

CALL HIST ( ITENP, A # B # NC , NRK , NP ) 

CALL HEANS( ITENP, XTEHP, NP, NRK, DDF ) 

CALL SUAP( -SP23.SU* , IER ) ; CALL LFF(XTENP,TENP2,TFNP3,NP,P,NR»: ) 
CALL SUAP ( ' SP24 . SU * , IER ) ; CALL HFF (XTENP, TENP2, TEhP3, NP , 0, NRK ) 
CALL DELETE C2RD') 

CALL RENAN ( -ZPP* , -ZRD* , IER ) 

CALL SIGNA(XTEHP,NRK,NP ) 

CALL AUTOCOR(XTENP,NP,NRK,N,RRH,P,Q > 

NP1-N+1 
DO 30 J-1.NP1 
RR(J)-RRH(J)+RR(J ) 

30 CONTINUE 

CALL PSDF(RRH,N,H,XL,GKC ) 

CALL INF<GKC,GKCP,G,H) 

CALL AC(GKCP,GKC3P -N,H,XA ) 

CALL SPDF(N,GKC3P '.H > 

DO 40 J-1,HP1 
GH( J)-GH(J)+GXH(J ) 

40 CONTINUE 

CALL DELETECRD* ) 


A-18 


20 


50 


100 

60 


70 

R 


CALL DELETECNRD* > 

CALI DELETE ( 'ZRO* ) 

CALL DELETECZP') 

CALL DELETECTD* > 

CALL F OPEN <3, *RD*,200O) 
call FOPEN(is,*sRD*.2ee0> 

IF ( I . GE . NR )GO TO 20 
Ii-I+1 

CALL READRC IS» II# ITEHP, 1 » IER ) 
CALL URITR(3,i» ITEfiP, 1 , IER ) 

CALL FCLOS CIS) 

CALL FCL0S(3 ) 

CONTINUE 

CALL DELETE (' PAR * ) 

DO 50 J-1,KP1 
RR< J )-RR( J )/NR 
GH< J )-GH( J )/NR 
CONTINUE 

TYPE* R AUERAGE R-HAT 


■ 


DO 60 J«1,KP1 


IR-J-1 

urite 

FORHAT (IX, 13, 8X, Ell. 4) 

CONTINUE 

TYPE* K AUERAGE GH(K)* 

DO 70 J ■ 1 , API 


IR-J-1 

URITEC 10, 100 ) IR.GH(J) 
CONTINUE 

CALL LAST (RR,GH,NP1 ,NR ) 


STOP 

END 


oonnooo 


TYPE SP26 

FILE SP26-LAST 
PROGRAMMER NAME t S.P. SINGH 
REUISIONi© 

DATE * 8/3/79 

PURPOSE * THIS MODULE COMPUTES STANDARD 
DEUIATION AND ERROR FOR THE AVERAGED 
R-HAT AND G-HAT ARRAYS. 

COMPILER NOSTACIC 
SUBROUTINE LAST ( RR. GH, HP1 , NR ) 

DIMENSION RR(101 >,GH( 101 ) 

RSQ-0. 

GSQ-0. 

DO 80 I-l.HPl 
RSQ-RSQMRR< I )**2 ) 

GSQ s GSO+GH( I )S>2 
80 CONTINUE 

RSD-fQRT(RSO/(MPl-l ) ) 
GSD-SGRT(GSQ/(MPl-i )) 

RN-NR 

ERH-RSDzSGRT(RN) 

EGH-GSD/SQRT(RN) 

TYPE 'STANDARD DEUIATION FOR R-HAT*, RSD 

TYPE 'STANDARD DEUIATION FOR G-HAT', GSD 

TYPE 'ERROR FOR R-HAT ',ERH 

TYPE 'ERROR FOR G-HAT', EGH 

STOP 

END 


4. SAMPLE RUNS 


Sample run results are included in the following 
pages for a file containing 3 records. Same data is used for 
each of the four cases namely without filters (SP1) , with low 
frequency filter only (SPL) , with high frequency filter only 
(SPH) and with both filters (SP2). 


RWRonn iHiuTY of mi': 

ORIGIN A I, I 1 \<; IS »’ \ 


FILE NUMRER,(**I „ 

DIFFUSE DENSITY FACTOR, DDF*1 

mxirum lac factor. R*« 

HEIGHT OF SCAfRUNC SLIT. XL- 1 
UIDTH OF SCAtRIINO SLIT.XA-l 
DATA PREPARATION COMMENCED 
FILE IS DE INC READ FRON TARE 
HUMIER OF RECORDS READ* _,_2 
FILE READ FRON TAPE AND WRITTEN 

ta W<D(k> .>o»0»ytS 

HOMIER OF SAMPLES- tRI 

X COORDINATE (NICRONS I- -3607 
V COORDINATE (HICROHS >• • 

DELTA X <M> IRICROMS** I* 

DATA IS IE I NO FORMATTED 
).M 3.3S }.)9 3.3i 3-35 

3 « 3.3* 322 3-» 3.37 

3.29 3.28 3.39 3S2 

3.4* 3.27 3.36 3.47 336 

3.42 3.43 3. 3D 3.29 3«» 

3.34 3.26 3.36 3. 2D 3. 3D 

3.31 3.37 3.43 3.(4 3.32 

3.31 3-41 3.37 3.27 3.47 

3.31 3.31 3.K 3. S3 3.43 

3.34 3.31 3.42 3.46 3.41 


ON DISK 


3.42 3. 3D 3-22 3. 29 33/ 3.*i ;»•« 

3. 2S 3.29 3.29 3.39 3.S2 3» 3.13 

3.44 3.27 3.36 3.47 3-36 3.36 3.50 

3 42 7 43 3 3D 3.29 3.41 3.33 3.33 

jS 3i3 3. !• 3 3* 3.24 3.34 

j!31 3.37 3.43 3.(4 3.32 3. 3D 3-39 

3.31 3.41 3.37 3.27 3.47 3.42 

3.31 3.31 3. 52 3. S3 3. 43 3.43 3.34 

3.34 3.31 3.42 3.46 3-41 3.29 3.46 

DATA* FORMATTING COMPUTED 
RAW DATA 6T0RED ON DISK. 

RECORD NO. 1 

HISTOGRAM DEINC COMPUTED 

TOTAL NO. OF CLASSES* IDD 

CLASS NO. HO. f r DATA POINTS 

79 * 

DD D 

«1 * 

12 3 

D3 22 

•4 23 

n id 

D6 16 

*7 4 

M 4 

■9 3 

RAW (CAM • D.26D5D6E 4 

REAM OF NORMALIZED DATA* D . 1DDDDDE t 
KJW SOUARE (MUJE • D.49T76DE -3 

ROOT hern SQUARE VALUE - 0.22 33 1 7E -• 

AUTOCORRELATION IS DEI NO COWED 

■ R-MAT NORN. R-MAT 


D.4D3DE -3 
D.113DE -3 
-0.6913E -4 
«u.«K> 
S.166CC -I 
D. 136 If -I 
D.4766C -2 
0KC3P/K 1 
D.1562E -1 
D.1353E -1 
D.4703E -2 


D.1DDDC 1 

D.22D9E D 
-D.14DDE D 
OKCJPir i 

D.IS62E -1 
0.13S3E -1 
0.4793E -3 
GH(K) 
D.14S7E -t 
D.UD7C -1 
0.9161E -2 
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A-26 


r 


mw newt • A.asaissc 4 

•CAN or NORMALIZED DATA* t.U MMI 1 
DATA I* KINO FILTERED FOR HIGH FREOUENCV 
ME AH MUARC VALUE • 4 

MOOT ICAN MUARC VALUE • -2 

AUTOCOMMELATIOH IS KINO CONFUTED 

R R-HAT NORM. R-HAT 

• t.13211 -3 •-!•••§ i 

A. HBOC -3 A.A96SE • 

2 A.97CSE -4 M.73SK • 


TIC KNOHINATOR OF SIlC FUNCTION KLOU 1.K-* 
K GXC<K> GM.CH*. > 

R A. 1 I2*E -I A.U24E -I 

1 R.R417E -3 A.2IA4E -I 

2 -A. 1D2IE -3 -A.4S52E -2 

«CF<K) G*C3PiK) 

A A.1124E -I A.I124E -I 

A.21A4E -I A.21A7E -I 

a -A.4552C -2 -A.457K -2 

OCC3R(Kl OH«CI 

A A.1124E -l A.IAIK -1 

A.21A7E -I A.122K -1 

2 -A.4STK -2 A.8247E -2 

R AUCRACE R-HAT 

A A.1491E -3 

| A.132K -3 

2 A.IA82C -3 

t AUCRACE CH(K> 

A A.I8S7C -1 

1 A. HUE -I 

a A.97A3E -a 

STANDARD DEVIATION FOR R-HAT A.1SA63AC -3 

STANDARD DEVIATION FOR C-HAT A.178743E -I 

ERROR FOR 8-HAT A.B27398E -4 

ERROR FOR C-HAT A.1A3I97E -1 

STOF 

■ 




FILE NUNKR.HF-1 
DIFFUSE WHS I TV FACTOR, DDF -1 
HAXIHUM LAO F ACTOR. A*2 
HEIGHT OF SCANNING SLIT.Xl-1 
WIDTH OF SCATWIM SUT.XA-1 
LOU FKOUENCV FILTER PARAKTCR.P-IO 
MICH FRCOUCNCV FILTER FAR ARE TER, 0*1 
DATA PREPARATION COHKHCED 
FILE IS KINO READ FRON TAPE 
NUHDER OF RECORDS READ* 3 

FILE READ FRON TAPE AND WRITTEN ON DISK 
Du klOik: ,>»>Osy >• 

HOMIER OF SAMPLES* ID! 

-3607 

D 

12 


-0.I17M -3 
QECP(K) 

-0.SD44C 
0KC3P(K ) 

-0 

0.S40K 

-a 

o.64aac 

■a 

0.3IE7E 

-l 

0.31S2E 

-i 

-0.S044C 

-a 

-O.SOTK 

-a 

CKC3P(K) 

GH(K > 

0.S462E 

-a 

O.ltSSC 

-i 

0.31S2C 

-i 

O.ISSK 

-i 

-O.SOTK 

-a 

0.I2S2E 

-i 


1 

« 

ECORD NO. 


HISTOOPNH KING COMPUTED 


> 

I 

to 

-J 


X COORDINATE (MICRONS )• 

V COORDINATE (MICRONS)* 
DELTA X (H) (MICRONS )• 
DATA IS KINO FORMATTED 
3.34 3.3S 3.30 3 “ 

3.42 3. 3D 3.22 3 

3.2S 3.20 3.20 3 
3.4# 3.2? 3.36 3 

3.42 3.43 3. 3D 3 
3.34 3.26 3.30 3 

3.31 3.3? 3.43 3 

3.31 3.40 3.37 3 

3.31 3.30 3.S2 3 

3.34 3.31 3.42 3 

3.40 

DATA FORMATTING COMPUTED 
RAU DATA STORED ON DISK. 

RECORD NO. 1 

HISTOGRAM KINS COMP U T E D 

TOTAL NO. OF CLASSES* 


CLASSES' 

NO. 


20 

30 

47 


14 

2? 


3? 

62 

36 

41 

30 

32 

47 

43 

41 


3.32 

3.41 
3.20 
3.36 

3.33 
3.24 
3.30 

3.42 

3.43 
3.20 


3.31 

3.20 

3.13 

3.S4 

3.33 

3.34 
3.30 
3.61 
3.34 
3.46 


30 3 

32 3 

22 3 
34 3 

43 3 
36 3 

34 3 

36 3 

33 3 
30 3 


35 

36 
34 
2? 
31 


34 
61 

35 

36 
31 
27 


41 


CLASS NO. 


too 


NO. or DATA POINTS 


TOTAL NO. OF 
CLASS NO. 

DO 
01 
•2 
03 
D4 

as 
06 
■7 
DS 
09 
9R 
91 

MU MEAN • 0.26D336C 4 

MEAN OF NORMALIZED DATA* 0.1 MM D C I 
MTA IS KINO FILTEKD FOR LOU FRCOUCNCV 
DAT* IS KINO FILTEKD FOR NIGH FKOUENCV 
MEAN SCUARE VALUE - 0.467603E 

HOOT KAN SOUAK VALUE • 0. 

AUTOCORRELATION IS KINO COMUTED 
R R-HAT NORN. 

0 0.1279C -3 0.14 

1 0.I944C -3 t.OISK 


OF MTA POINTS 
1 
0 
10 
13 
10 
22 
12 
D 
3 

a 

i 
t 


-2 


70 

a 

2 

0.7I63C -4 

• o 

M 

0 

THE 

DENOMINATOR OF SIK 

FUNCTION iCLOU 1 .PC -06 

01 

a 

C 

CKC(K) 

2(1) 

t2 

7 

0 

O.08MC -2 

0-98MC -2 

03 

32 

1 

D. 1 3521 -2 

D.3379C -1 

■4 

23 

2 

-0.221 II -3 

-O.S52K -2 

K 

ia 

E 

GKCP(K) 

CEC3P(K) 

>6 

16 

0 

0.9BDK -2 

0.988K -2 

■7 

4 

1 

P.3379C -1 

0.3384C -1 

M 

4 

2 

-0.6S2K -2 

-O.S56K -2 


ENCOUNTERED IN IK 


MU MEAN • 

KAN OT NORMALIZED DATA* 0.11 
MTA 10 KINO FILTERED FOR LOU FKOUENCV 
MTA 10 KIM FILTERED FOR HIGH FKOUENCV 
KAN SOUAK VALUE • 0.SI7979C -4 

ROOT HERN S OUA K VALUE * 0.7004621 -2 

AUTO C O RK LATION IS KINO COMPUTED 

■ R-HAT NORM. R-HAT 

0 0.2994C -4 9. tOOCC 1 

1 9.S96K -4 0./361C 0 

■ D.2D32E -4 0.3499C 0 

TM DENOMINATOR OF OIK FUNCTION KLOU t.K-M ENCOUNTERED IN IMF 


6EC3P(K I 

0.9000E -2 

0.33D4C -1 

-2 


K 

0 

1 

2 


KCORD NO. 3 

HISTOGRAM KING COMPUTED 
TOTAL NO. OF CLASSES* 


GM(K> 
0.2102C -1 
0.1 TORE -1 
0.1414C -1 


K 

D 

1 


GKC(K) 

0.640K -a 

0.1263C -a 


GKCP(K) 
0.6462E -2 
0.3167C -I 


CLASS NO. 

70 

DO 
tt 

•a 

S3 


NO. OF DATA POINTS 


* 

m 



A-28 


TYPE SP16 

C FILE SP16-AVG 

C PROGRAMMER NAME: S.P. SINGH 

C REVISION i 0 

C DATE* 3/S/79 

C PURPOSE: THIS MODULE COMPUTES THE RAU MEAN OF THE DATA. 

COMPILER NOSTACIC 

SUBROUTINE AVG< ITEHP, ICHR,NR,NP,RAUM ) 

C SP1S IS AUG 

REAL ISUM 

DIMENSION ITEMP( 1606 ) # NP(50 ) 

CALL FOPENC 1CHR, “RD* ,2060 ) 

INP-0 

ISUM-0 

C 

C SUM ALL THE POINTS IN THE FILE 

C 

DO 10 I « 1 , NR 

CALL READRilCHR.I, I TEMP, 1 , IER ) 

NPI -NP( I ) 

DO 20 J-1,NPI 
ISUH-ISUH+ITENPU) 

20 CONTINUE 

INF-INP+NPI 
10 CONTINUE 

RAUM-ISUM/INP 

TYPE *RAU MEAN - *,RAUM 

CALL FCLOS(ICHR) 

RETURN 

END 

R 


